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In  support  of  the  Defense  Nuclear  Agency’s  HYDRO  PLUS  program,  SRI  performed  a  series  of  uniaxial  strain 
experiments  using  in-contact  explosive  techniques  to  investigate  the  effects  of  in  situ  and  artificial  joints  on  wave 
propagation  in  calcium  carbonate  (CaCOa)  rock.  The  in  situ  Joints  investigated  were  healed,  calcite-filled  fissures 
in  a  nonporous  limestone  matrix.  These  joints  did  not  have  a  significant  effect  on  the  observed  response.  On  the 
other  hand,  artificial,  sand-filled  joints  in  marble  caused  an  increase  in  the  peak  particle  velocity  upstream  from  the 
joint  and  a  delayed  time  of  arrival  (TOA)  of  the  wave  at  the  gage  downstream  from  the  joint.  Both  peak  particle 
velocity  and  TOA  are  used  within  the  framework  of  the  HYDRO  PLUS  yield  estimation  methodology. 

A  spherical  wave  experiment  was  also  conducted  on  a  sample  of  intact  Danby  marble  to  complement  the  uniaxial 
strain  experiments  and  to  provide  data  that  can  be  used  to  validate  and/or  calibrate  constitutive  models  for  use  in 
finite  element  and  finite  difference  hydrocodes.  The  results  of  this  experiment  indicate  that  the  occurrence  of 
cracking  can  have  a  significant  effect  on  the  observed  response  if  measurements  are  made  in  the  strength- 
dominated  region. 
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EXECUTIVE  SUMMARY 


The  present  investigation  was  undertaken  as  part  of  the  HYDRO  PLUS  program  to  gain  an  improved 
knowledge  of  the  material  behavior  parameters  and  thus  increase  confidence  in  the  HYDRO  PLUS  yield 
estimation  methodology.  For  this  purpose,  a  series  of  uniaxial  strain  experiments  were  performed  using 
in-contact  explosive  techniques  to  investigate  the  effect  of  in  situ  and  artificial  joints  on  wave 
propagation  in  calcium  carbonate  (CaC03)  rock.  A  spherical  wave  experiment  was  also  conducted  on  a 
sample  of  intact  Danby  marble  to  complement  the  uniaxial  strain  experiments  and  provide  data  that  can 
be  used  to  validate  and/or  calibrate  constitutive  models  for  use  in  tinite  element  and  finite  difference 
hydrocodes.  Hydrocode  simulations  were  used  to  aid  in  the  interpretation  of  the  experimental  data.  The 
uniaxial  strain  experiments  were  simulated  using  the  one-dimensional  finite  difference  code  SRI  PUFF, 
and  the  spherical  wave  experiment  was  simulated  in  the  finite  element  code  DYNA2D  using  a 
multiplane  cracking  model. 

The  in  situ  joints  investigated  were  healed,  calcite-filled  fissures  in  a  limestone  matrix,  Utah  Test  and 
Training  Range  (UTTR)  limestone,  composed  primarily  of  calcite  (62%)  and  quartz  (31%)  with  traces  of 
dolomite  (2%),  K-feldspar  (2%),  illite  (2%),  and  pyrite  (1%).  These  in  situ  joints  did  not  have  a 
significant  effect  on  the  observed  response. 

The  artificial  joints  were  embedded  in  Danby  marble  specimens.  The  joint  parameters  varied  during  our 
investigation  included  surface  texture  (smooth  versus  rough),  thickness  (closed  versus  sand-filled), 
number  (single  versus  triple),  and  orientation  (normal  versus  inclined  with  respect  to  the  direction  of 
shock  propagation).  Intact  marble  specimens  were  also  tested  to  provide  baseline  data  for  comparison 
with  data  from  the  jointed  samples. 

Closed  joints  did  not  have  a  measurable  effect  on  the  behavior  of  the  rock,  nor  did  the  surface  texture  of 
the  joint.  However,  significant  effects  were  produced  when  a  sand-filled  joint  was  introduced  into  the 
sample.  Upstream  from  the  joint,  a  temporary  increase  in  the  magnitude  of  the  particle  velocity  profiles 
was  observed.  Downstream  from  the  joint,  the  effect  of  the  joint  was  manifested  as  a  delayed  time  of 
arrival  (TOA)  of  the  wave  at  the  gage  plane.  The  perturbations  to  the  free  field  were  greatest  for  the 
triple  normal  joint  configuration  and  smaller  for  the  single  inclined  joint  configuration,  while  the  single 
normal  joint  configuration  had  the  smallest  measurable  effect. 

From  these  results,  we  conclude  that  the  presence  of  sand-filled  joints  can  have  a  significant  effect  on 
the  HYDRO  PLUS  yield  estimation  methodology,  because  sand-filled  joints  perturb  the  free-field 
motion  in  a  manner  that  causes  an  increase  in  the  peak  particle  velocity  upstream  from  the  joint  and  a 
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delayed  TOA  downstream  from  the  joint.  Both  peak  particle  velocity  and  TOA  are  used  within  the 
framework  of  the  HYDRO  PLUS  yield  estimation  methodology. 

In  addition  to  the  effects  of  joints  on  wave  propagation  described  above,  the  results  of  the  uniaxial  strain 
experiments  yielded  information  regarding  the  precursor  and  a  high  pressure  phase  transformation  in 
CaCOs  rock.  The  precursor  was  observed  both  during  loading  and  unloading  and  was  most  likely 
associated  with  the  I-II  and  H-III  phase  ttansitions  in  calcite.  Based  on  the  loading  data,  the  magnitude 
of  the  precursor  was  estimated  to  be  between  16  and  19  kbar  in  marble  and  14.3  kbar  in  UTTR 
limestone.  The  high  pressure  phase  transition  was  observed  in  the  marble  experiments.  The  stress  level 
associated  with  this  transition  was  about  118  kbar. 

In  the  spherical  wave  experiment  on  Danby  marble,  extensive  radial  and  circumferential  crack  networks 
were  observed  in  the  specimen  after  the  test.  The  circumferential  crack  patterns  were  asymmetric  and 
led  to  load-induced  anisotropy  in  the  late-time  response  of  the  sample,  as  indicated  by  the  systematic 
differences  between  the  particle  velocity  histories  recorded  on  diametrically  opposite  sides  of  the 
explosive  charge.  Since  the  initial  loading  was  spherically  symmetric,  the  lack  of  symmetry  in  the  late¬ 
time  particle  velocity  histories  was  attributed  to  asymmetric  cracking.  Code  simulation  results  indicated 
that  close  to  the  charge  cavity  the  behavior  of  the  rock  was  dominated  by  porous  compaction,  whereas 
away  from  the  charge  the  behavior  was  dominated  by  crack  propagation. 

Accurate  modeling  of  the  cracking  process  thus  appears  to  be  necessary  to  allow  assessment  of  the 
potential  impact  of  cracking  on  the  HYDRO  PLUS  yield  estimation  methodology  if  measurements  are 
made  in  the  strength-dominated  region. 
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♦The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  =  1  event/s. 
**The  Gray  (Gy)  is  fiie  SI  unit  of  absorbed  radiation. 
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SECTION  1 
OVERVIEW 


1.1  BACKGROUND  AND  OBJECTIVES. 

The  present  investigation  was  undertaken  as  part  of  the  HYDRO  PLUS  program  to  gain  an  improved 
knowledge  of  the  material  behavior  parameters  and  thus  increase  confidence  in  the  HYDRO  PLUS  yield 
estimation  methodology.  With  this  methodology,  the  yield  of  a  large,  nonstandard,  underground  nuclear 
event  is  estimated  based  on  a  combination  of  experimental  measurements  of  the  shock  wave  parameters 
(stress,  particle  velocity,  shock  velocity)  near  the  source  and  computational  predictions  of  the  close-in 
ground  motion  as  a  function  of  the  yield  of  the  nuclear  event.  The  reliability  of  the  estimated  yield  is 
therefore  closely  related  to  the  quality  of  the  experimental  ground  shock  data  and  the  accuracy  of  the 
constitutive  model  used  in  the  numerical  simulations. 

During  the  early  stages  of  the  investigation,  SRI  performed  a  literature  review  to  identify  regions  of  the 
high  pressure  equation  of  state  of  rock  that  are  of  interest  to  the  HYDRO  PLUS  program  and  have  not 
been  explored  experimentally.  The  material  behavior  parameters  emphasized  in  this  review  were  the 
loading  Hugoniot,  the  unloading  adiabat,  phase  transitions,  and  the  effect  of  joints  on  wave  propagation. 
Our  findings,  which  are  detailed  in  a  later  section,  revealed  that  the  effect  of  joints  on  wave  propagation 
is  the  least  investigated  aspect  of  the  behavior  of  rock.  For  this  reason,  investigating  the  effect  of  joints 
on  wave  propagation  became  the  primary  objective  of  the  present  investigation.  Calcite  (calcium 
carbonate,  chemical  formula  CaCOs)  was  chosen  for  this  study;  partly  because  it  is  one  of  the  most 
abundant  carbonate  mineral  types  in  the  crust  of  the  earth  and  partly  because  of  DNA's  ongoing  interest 
in  obtaining  statistical  variations  in  material  properties  for  jointed  limestone  (a  calcite  rxx:k)  and  in 
developing  techniques  for  including  the  observed  effects  in  material  models  for  ground  shock. 

The  experiments  performed  during  the  study  are  summarized  in  Table  1-1.  As  indicated,  all  but  one  of 
the  experiments  were  performed  under  uniaxial  strain  loading  conditions.  Three  of  the  uniaxial  strain 
experiments  were  performed  on  UTTR  limestone  samples  with  naturally  occurring,  in  situ  joints.  Danby 
(Vermont)  marble  was  used  in  all  the  remaining  experiments.  Well  characterized  artificial  joints  of 
various  configurations  were  introduced  into  the  marble  specimens  to  investigate  the  effect  of  joints  on 
wave  propagation.  A  spherical  wave  experiment  was  also  conducted  on  a  sample  of  intact  Danby 
marble  to  complement  the  uniaxial  strain  experiments  and  to  provide  data  that  can  be  used  to  validate 
and/or  calibrate  constitutive  models  for  use  in  finite  element  and  finite  difference  hydrxjcodes. 
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1.2  SUMMARY. 


1.2.1  Literature  Review. 

Numerous  studies  have  been  conducted  over  the  past  30  years  to  characterize  the  behavior  of  calcium 
carbonate  rock  under  dynamic  loadings.  Many  of  the  early  dynamic  wave  propagation  studies 
emphasized  the  characteristics  of  the  loading  Hugoniot  and,  to  a  lesser  extent,  the  unloading  adiabat  and 
the  effects  of  strength,  porosity,  and  phase  transformations  on  the  material  behavior. 

Static  high  pressure  studies  played  a  significant  role  in  furthering  our  understanding  of  the  behavior  of 
calcium  carbonate  rock.  The  pressure  and  volume  change  associated  with  the  calcite  I  -  calcite  11  and 
calcite  11  -  calcite  IE  phase  transitions  were  first  measured  accurately  in  static  experiments.  These 
transitions  have  since  been  observed  in  shock  wave  propagation  experiments  and  are  now  known  to 
affect  the  shape  of  both  the  loading  Hugoniot  and  the  unloading  adiabat  of  calcite  rock  in  a  predictable 
fashion.  The  loading  and  unloading  behaviors  of  calcite  rock  are  also  affected  by  another  phase 
transformation  -  the  calcite  HI  -  calcite  IV  transition.  This  transition  occurs  at  stresses  in  excess  of 
100  kbar  and  has  only  been  observed  in  dynamic  wave  propagation  experiments.  Summaries  of  static 
and  dynamic  phase  transformation  measurements  in  calcium  carbonate  rock  are  provided  in  Tables  1-2 
and  1-3,  respectively. 

The  Hugoniot  of  calcium  carbonate  rock  has  been  investigated  for  stresses  up  to  about  500  kbar.  For 
nonporous  rock,  the  Hugoniot  is  well  characterized  and  results  from  numerous  investigations  are  very 
reproducible,  as  depicted  in  Figure  1-1.  The  Hugoniot  elastic  limit  (HEL)  of  calcite  occurs  in  the  stress 
range  from  10  to  20  kbar  and  depends  on  the  porosity  and  graininess  of  the  sample.  The  HEL  is 
associated  with  the  first  two  polymorphic  calcite  transitions  (I  -  II  and  11  -  HI). 

The  unloading  adiabat  is  not  too  far  removed  from  the  loading  Hugoniot,  partly  because  of  the  fully 
reversible  nature  of  the  calcite  phase  transitions.  In  nonporous  calcite,  little  hysteresis  is  observed  upon 
complete  unloading  from  a  shock-loaded  state.  As  one  might  expect,  hysteresis  increases  with 
increasing  porosity. 

The  effect  of  joints  on  material  response  is  the  least  investigated  aspect  of  the  behavior  of  calcite  rock. 
Our  literature  review  unveiled  only  two  studies  aimed  at  investigating  the  effect  of  joints  on  wave 
propagation.  These  studies  showed  drastic  differences  between  the  responses  of  intact  and  jointed 
samples.  This  discovery,  combined  with  DNA's  interest  in  obtaining  statistical  variations  in  material 
properties  for  jointed  limestone  and  developing  techniques  for  including  the  observed  effects  in  material 
models  for  ground  shock,  led  to  the  consensus  that  the  experiments  to  be  conducted  in  this  study  should 
emphasize  the  effect  of  joints  on  material  behavior. 
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Table  1-3.  Dynamic  experimental  measurements  of  the  pressures  associated 
with  the  various  calcite  phase  transitions. 
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Figure  1-1 .  A  collection  of  Hugoniot  data  for  nonporous  (less  than  0.5%  porosity)  calcium  carbonate  (CaCOg)  rock 
in  the  stress-volumetric  strain  plane. 


1.2.2  Uniaxial  Strain  Experiments. 

A  series  of  uniaxial  strain  experiments  was  performed  using  in-contact  explosive  techniques  to 
investigate  the  effect  of  in  situ  and  artiricial  joints  on  wave  propagation  in  calcium  carbonate  (CaC03) 
rock.  The  experiments  were  motivated  by  a  lack  of  experimental  data  on  wave  propagation  in  jointed 
rock  and  were  aimed  at  providing  results  that  will  further  our  understanding  of  the  effect  of  joints  on 
wave  propagation  and  provide  a  consistent  set  of  results  that  can  be  used  to  develop,  calibrate,  or 
validate  constitutive  models  suitable  for  use  in  hydrocode  calculations. 

The  study  consisted  of  experimental  and  analytical  investigations.  Experimentally,  uniaxial  strain 
experiments  were  conducted  on  both  intact  and  jointed  calcite  rock  samples.  The  experimental 
configurations  investigated  are  shown  in  Figure  1-2.  In  the  experiments,  in-contact  explosives  initiated 
using  15.24-cm-diameter  (6.0  in.)  plane  wave  lenses  were  used  to  provide  shock  loadings  in  the  rock 
samples.  Several  specimen  designs  incorporating  different  joint  configurations  were  investigated. 
Experimental  diagnostics  consisted  of  in-material  Lagrangian  particle  velocity  gages  embedded  between 
rock  layers  at  several  serial  locations  within  the  specimen.  Duplicate  particle  velocity  measurements 
were  obtained  on  each  of  four  gage  planes.  The  particle  velocity  gages  on  the  first  gage  plane  (nearest 
to  the  explosives)  were  used  to  record  the  oncoming  stress  wave  generated  by  the  explosives  and  thereby 
provide  a  basis  for  comparing  the  results  of  different  experiments  with  different  joint  configurations. 

The  remaining  gages  recorded  the  wave  profile  after  passage  of  the  wave  through  the  jointed  section  of 
the  specimen. 

Analytically,  Lagrangian  finite  difference  hydrocode  simulations  were  carried  out  to  help  in  the 
interpretation  of  experimental  results  and  provide  added  insight  into  the  effect  of  joints  on  the  response 
of  calcium  carbonate  rock. 

1.2.2.1  Effect  of  .Joints  on  Wave  Propagation.  Two  types  of  joints  were  investigated:  in  situ  joints  in 
UTTR  limestone,  and  artificial  joints  in  Danby  (Vermont)  marble.  The  in  situ  joints  investigated  were 
healed,  calcite-filled  fissures  in  a  UTTR  limestone  matrix  composed  of  calcite  (62%)  and  quartz  (31%) 
with  traces  of  dolomite  (2%),  K-feldspar  (2%),  illite  (2%),  and  pyrite  (1%).  Our  limited  results  from  the 
UTTR  limestone  experiments  indicated  that  the  in  situ  joints  did  not  have  a  significant  effect  on  the 
observed  response. 

The  artificial  joints  were  embedded  in  Danby  marble  specimens.  The  joint  parameters  that  were  varied 
included  surface  texture  (smooth  versus  rough),  thickness  (closed  versus  2-mm-thick  sand-filled), 
number  (single  versus  triple),  and  orientation  (normal  versus  inclined  with  respect  to  the  direction  of 
shock  propagation).  Intact  marble  specimens  were  also  tested  to  provide  baseline  data  for  comparison 
with  data  for  the  jointed  sample. 
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(a)  Configuration  for  experiments  on  intact  rock  (basic  configuration). 


(b)  Modification  of  the  basic  configuration  to  include  a  single,  normal  joint. 


(c)  Modification  of  the  basic  configuration  to  include  a  triple,  normal  joint. 


(d)  Modification  of  the  basic  configuration  to  include  a  single,  inclined  joint. 


Figure  1-2.  Configurations  of  the  uniaxial  strain  experiments. 


Closed  joints  did  not  have  a  measurable  effect  on  the  behavior  of  the  rock,  nor  did  the  surface  texture  of 
the  joint.  However,  significant  effects  were  produced  when  a  sand-tilled  joint  was  introduced  into  the 
sample.  The  perturbations  to  the  free  field  were  greatest  for  the  triple  normal  joint  configuration  and 
smaller  for  the  single  inclined  joint  configuration,  while  the  single  normal  joint  configuration  had  the 
smallest  measurable  effect. 

To  better  understand  the  effect  of  artificial  joints  on  wave  propagation,  we  performed  a  series  of 
hydrocode  simulations  using  the  finite  difference,  wave  propagation  hydrocode  SRI  PUFF  (Seaman  and 
Curran,  1978).  First,  we  simulated  the  response  of  an  intact  marble  specimen  having  the  configuration 
shown  in  Figure  l-2(a).  We  then  simulated  the  response  of  jointed  specimens  having  the  configurations 
shown  in  Figures  l-2(b)  (single  joint)  and  l-2(c)  (triple  joint).  The  effect  of  joints  was  established  by 
comparing  the  results  of  the  benchmark  simulation  in  intact  marble  with  the  results  of  simulations  in 
jointed  marble. 

Figure  1-3  shows  the  simulated  particle  velocity  histories  at  each  of  the  gage  planes  and  at  the  fiee 
surface  of  the  intact  marble  experiment.  The  first  peak  in  the  particle  velocity  histories  corresponds  to 
peak  stress,  and  it  is  the  same  on  all  the  gage  planes.  This  peak,  which  also  corresponds  to  the 
Chapman-Jouguet  (C-J)  pressure  of  the  explosives,  is  followed  by  a  gentle  unloading  wave  that  causes  a 
gradual  decrease  in  the  particle  velocity  profiles.  This  is  the  Taylor  unloading  wave,  which  originates  in 
the  explosive  gases.  The  second  plateau  in  the  particle  velocity  records  is  associated  with  the  stress-free 
state  in  the  specimen.  It  is  earned  through  the  specimen  by  a  backward-travelling  release  wave  that 
originates  at  the  back  surface  of  the  specimen. 

To  gain  insight  into  the  effect  of  joints  on  wave  propagation,  the  two  experimental  contiguradons  with 
normally  oriented  artificial  joints,  shown  in  Figure  l-2(b)  and  l-2(c),  were  also  simulated.  The  resulting 
particle  velocity  histories  at  the  first  and  second  gage  planes  are  compared  with  the  results  of  the  intact 
marble  simulation  in  Figures  1-4  and  1-5,  respectively.  Two  features  associated  with  the  effect  of  joints 
on  the  simulated  particle  velocity  profiles  are  readily  apparent.  The  first  of  these  features  is  depicted  in 
Figure  1-4:  The  particle  velocity  profiles  on  the  first  gage  plane  (ahead  of  the  joint)  are  modified  by 
spikes — one  spike  for  the  single  joint  configuration  and  three  spikes  for  the  triple  joint  configuration. 
These  spikes,  and  the  corresponding  decrease  in  stress,  result  from  wave  reflections  at  the  rock/joint 
interfaces.  Interestingly,  these  spikes  are  not  quite  as  apparent  on  the  gage  planes  located  behind  the 
joint,  although  on  the  second  gage  plane  the  velocity  history  associated  with  the  triple  joint 
configuration  shows  characteristics  similar  to  those  observed  on  the  first  gage  plane.  The  second  feature 
associated  with  the  effect  of  joints  on  the  simulated  particle  velocity  profiles  is  the  delayed  time  of 
arrival  (TOA)  of  the  wave  at  the  gage  planes  behind  the  joint.  This  delayed  TOA  is  barely  perceptible 
for  the  single  joint  configuration  but  more  pronounced  for  the  triple  joint  configuration.  The  TOA  delay 
can  be  related  to  the  number  of  reverberations  required  to  bring  the  joint  to  pressure  and  particle  velocity 
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PARTICLE  VELOCITY  (mm/ns) 


Figure  1  -3.  Simulated  particle  velocity  histories  for  a  typical  experiment  on 
intact  marble. 
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equilibrium  with  the  test  medium.  It  is  not  surprising,  therefore,  that  the  TOA  delay  associated  with  the 
triple  joint  is  approximately  three  times  that  associated  with  the  single  joint. 

The  same  effects  observed  in  the  PUFF  simulations  can  also  be  observed  in  the  experimental  data,  as 
illustrated  in  Figures  1-6  and  1-7  for  the  single  joint  and  triple  joint  configurations,  respectively.  In  each 
figure,  the  results  from  an  experiment  on  intact  marble  are  superimposed  onto  the  results  of  an 
experiment  on  jointed  marble.  Part  (a)  in  each  figure  corresponds  to  the  gage  location  in  front  of  the 
jointed  section  of  the  specimen,  whereas  part  (b)  corresponds  to  the  gage  location  behind  the  jointed 
section.  The  narrow,  large  amplitude  spikes  at  t  =  8  |is  in  Figure  l-6(a)  and  at  t  =  1 1  ^is  in  Figure  l-7(a) 
are  most  likely  due  to  an  air  gap  at  the  interface  between  the  joint  and  the  rock.  This  explanation  is 
plausible  because  the  spike  is  not  observed  at  any  of  the  gage  planes  behind  the  joint;  and  also  because 
the  spike  reaches  the  gage  before  the  expected  reflection  from  the  joint,  a  result  indicating  that  the 
source  of  the  spike  is  located  at  the  interface  between  the  gage  plane  and  the  joint. 

1.2.2.2  Implications  of  the  Results  on  the  HYDRO  PLUS  Program.  The  HYDRO  PLUS 
methodology  estimates  the  yield  of  a  nuclear  event  on  the  basis  of  a  combination  of  experimental 
measurements  of  the  shock  wave  parameters  near  the  source  and  computational  predictions  of  the  close- 
in  ground  motion.  The  methodology  relies  on  measurements  of  peak  stress,  peak  particle  velocity,  and 
shock  TOA  (and  hence  shock  velocity)  at  different  radii  from  the  charge.  The  presence  of  joints  was 
observed  to  have  a  wide  range  of  effects  on  the  magnitude  of  the  particle  velocity  and  the  TOA 
measured  in  our  uniaxial  strain  experiments.  On  one  hand,  negligible  effects  were  observed  when  the 
joint  was  either  closed  or  filled  with  a  material  that  has  an  impedance  similar  to  that  of  the  host  medium 
(i.e.,  calcite-filled  joints  in  UTTR  limestone).  On  the  other  hand,  significant  effects  were  observed  when 
the  joint  was  filled  with  a  material  that  has  an  impedance  substantially  different  from  that  of  the  host 
medium  (i.e.,  sand-filled  joints  in  marble). 

Figures  1-4, 1-6,  and  1-7  showed  that  a  sand-filled  joint  causes  a  temporary  increase  in  the  magnitude  of 
the  particle  velocity  upstream  from  the  joint.  The  TOA  of  the  joint  signature  at  the  gage  plane  depends 
on  the  distance  separating  the  joint  from  the  gage.  If  this  distance  is  small,  the  joint  signature  could 
affect  the  peak  of  the  particle  velocity  profile,  which  in  turn  may  lead  to  an  erroneous  HYDRO  PLUS 
yield  estimate.  It  is  possible  to  recognize  the  signature  of  a  joint  by  a  careful,  systematic  analysis  of  the 
recorded  data  (DeCarli,  1993),  but  it  is  not  easy  to  determine  the  free  field  velocity  from  the  joint- 
affected,  measured  velocity.  This  difficulty  may  have  an  adverse  effect  on  the  overall  precision  of  the 
yield  estimated  using  the  HYDRO  PLUS  methodology. 

Downstream  from  the  joint,  the  only  significant  effect  observed  in  our  uniaxial  strain  experiments  was 
delayed  TOA  of  the  shock  wave  at  the  gage  plane.  The  measured  peak  particle  velocity  (and  the 
simulated  peak  stress)  was  not  influenced  by  the  joint.  Thus,  the  presence  of  a  joint  upstream  from  the 
gage  would  only  have  a  significant  effect  on  the  HYDRO  PLUS  yield  estimation  methodology  if  the 
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(a)  Particle  velocity  histories  at  the  first  gage  plane. 


(b)  Particle  velocity  histories  at  the  second  gage  plane. 


Figure  1  -6.  Effect  of  a  single  sand-filled  joint  on  wave 
propagation. 
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(a)  Particle  velocity  histories  at  the  first  gage  plane. 


(b)  Particle  velocity  histories  at  the  second  gage  plane. 


Figure  1  -7.  Effect  of  a  triple  sand-filled  joint  on  wave 
propagation. 
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shock  velocity  calculated  based  on  TOA  measurements  is  used  to  estimate  the  yield.  It  is  worth  noting 
that  this  effect  becomes  more  significant  with  decreasing  distance  between  TOA  recording  stations. 

1.2.23  Hueoniot  and  Phase  Transitions.  In  addition  to  the  effects  of  joints  on  wave  propagation 
described  above,  the  results  of  the  uniaxial  strain  experiments  yielded  information  regarding  the 
precursor  and  a  high  pressure  phase  transformation  in  CaC03  rock.  The  precursor  was  observed  during 
both  loading  and  unloading  and  was  most  likely  associated  with  the  HI  and  Il-in  phase  transitions  in 
calcite.  Based  on  the  loading  data,  the  magnitude  of  the  precursor  was  estimated  to  be  between  16  and 
19  kbar  in  marble  and  14.3  kbar  in  UTTR  limestone.  The  high  pressure  phase  transition  was  observed  in 
the  marble  experiments.  The  stress  level  associated  with  this  transition  was  about  118  kbar. 

1.2.3  Spherical  Wave  Experiment. 

A  spherical  wave  experiment  was  conducted  on  a  sample  of  intact  Danby  marble  to  complement  the 
uniaxial  strain  experiments  and  provide  data  that  can  be  used  to  validate  and/or  calibrate  constitutive 
models  for  use  in  finite  element  and  finite  difference  hydrocodes.  In  the  experiment,  a  2-cm-diameter, 
6-g  EL-506D  (Detasheet)  explosive  charge  was  detonated  at  the  center  of  a  27-cm-diameter,  27-cm-long 
cylindrical  block  of  marble  as  shown  in  Figure  1-8.  The  specimen  assembly  was  placed  inside  a 
pressure  vessel  and  surrounded  by  a  solenoid  that  applied  a  constant,  vertically  oriented  magnetic  field 
through  the  particle  velocity  gage  planes.  A  6.9  MPa  static  overburden  pressure  was  applied  to  the 
specimen  during  the  experiment. 

Particle  velocity  histories  at  several  radii  from  the  charge  were  measured  by  monitoring  the  motion  of 
embedded  wire  gages  through  the  externally  applied  magnetic  field.  The  particle  velocity  gages 
consisted  of  strips  of  copper  wire  mounted  into  concentric,  precision-machined  grooves  in  the  specimen. 
As  shown  in  Figure  1-9,  two  diametrically  opposite  gages  were  used  at  each  of  eight  radii  starting  at  the 
charge  boundary  and  extending  outward  to  a  radius  of  5.0  cm.  In  most  applications,  the  active  element 
of  each  gage  would  cover  the  full  circumference  of  one  of  the  concentric  circles  shown  in  the  figure. 
However,  in  the  present  study  a  different  approach  was  taken  to  ensure  that  the  direction  along  which  the 
particle  velocity  histories  were  measured  coincided  with  the  direction  of  wave  propagation  in  the 
uniaxial  strain  experiments.  This  arrangement  allows  for  a  direct  comparison  between  the  results  of  the 
uniaxial  strain  experiments  and  the  results  of  the  spherical  wave  experiment  even  if  the  rock  behaves  in 
an  anisotropic  fashion. 

The  particle  velocity  histories  measured  in  the  experiment  on  Danby  marble  are  shown  in  Figure  1-10. 
Each  of  Figures  1- 10(a)  through  1- 10(h)  shows  the  two  particle  velocity  histories  measured  at  a  given 
range.  The  letters  ‘a’  and  ‘b’  (following  the  gage  number)  in  Figure  1-10  designate  the  two  gage  arrays 
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Figure  1-8.  Overall  configuration  of  the  spherical  wave  experiment  in  Danby  marble. 
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Figure  1-9.  Layout  of  the  particie  veiocity  gages  in  the  sphericai  wave  experiments. 
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on  diametrically  opposite  sides  of  the  charge.  The  particle  velocity  histories  at  the  two  ranges  closest  to 
the  charge  show  characteristics  that  are  different  from  the  remainder  of  the  data.  We  believe  this  is  the 
result  of  premature  gage  failure  due  to  the  high  stresses  generated  by  the  explosive  charge. 

The  particle  velocity  gages  at  the  2.0-cm  range  [Figure  l-10(c)]  captured  the  early- time  response  of  the 
rock  but  failed  to  record  the  subsequent  rebound  signal.  At  ranges  of  2.5  cm  or  larger,  the  particle 
velocity  histories  in  the  specimen  were  recorded  for  the  duration  of  the  experiment.  At  early  times, 
these  histories  are  characterized  by  a  shaip  rise  to  peak  followed  by  an  outward  motion  that  lasts  for 
about  10  \is.  This  early-time  response  of  the  rock  is  very  reproducible,  as  indicated  by  the  nearly 

identical  responses  of  the  gage  arrays  on  either  side  of  the  explosive  charge  [Figures  l-lO(c)  through 
l-10(h)]. 

The  outward  expansion  phase  is  followed  by  a  rebound  phase.  In  previous  experiments  on  nonporous 
rock  (Miller  and  Florence,  1990)  and  on  porous  CaCOs  rock  (Gefken  and  Rorence,  1992),  the  rebound 
signal  caused  an  inward  contraction  in  the  motion  of  the  rock  sample  that  was  manifested  by  a  reversal 
of  the  sign  of  the  particle  velocity  histories.  In  the  present  experiment,  the  particle  velocities  recorded  at 
all  ranges  remained  positive  throughout  the  rebound  phase  of  the  signal.  This  difference  between  the 
rebound  signals  measured  in  the  present  experiment  and  those  measured  in  other  investigations  can  be 
related  to  cracking  of  the  rock  sample.  Extensive  cracking  was  not  observed  in  either  of  the  two  studies 
cited  above.  In  contrast,  the  Danby  marble  specimen  of  the  present  study  was  severely  cracked,  as 
shown  in  Figure  1-11.  Two  distinct  networks  of  cracks  can  be  seen  in  each  half  of  the  specimen.  The 
first  network  consists  of  numerous  cracks  emanating  from  the  charge  boundary  and  extending  radially 
outward  toward  the  free  surface  of  the  specimen.  These  radial  cracks  extend  throughout  the  gaged 
region  of  the  specimen,  and  a  few  cracks  propagate  all  the  way  to  the  free  surface.  The  second  network 
of  cracks  consists  primarily  of  circumferential  cracks  and  does  not  appear  to  be  symmetric  with  respect 
to  the  center  of  the  explosive  charge.  An  additional  crack  that  does  not  appear  to  be  associated  with 
either  of  the  crack  networks  mentioned  so  far  can  be  observed  spanning  the  whole  specimen  surface. 

This  crack  follows  the  path  of  a  pre-existing  in  5/m  joint  in  the  rock. 

We  believe  that  the  extensive  cracking  in  the  Danby  marble  specimen  caused  degradation  of  the  material 
properties  and  reduced  the  stiffness  of  the  rock.  For  this  reason,  the  specimen  unloaded  along  a 
“damaged”  modulus  that  was  substantially  lower  than  the  intact  modulus.  Consequently,  much  of  the 
early-time  deformation  was  not  recovered  during  the  contraction,  or  rebound,  phase. 

The  results  of  the  spherical  wave  experiment  also  show  evidence  of  load-induced  anisotropy,  a  response 
mode  usually  associated  with  cracking  in  brittle  materials.  This  anisotropy  is  manifested  as  a  difference 
in  the  recorded  particle  velocity  histories  along  the  two  different  gage  arrays  of  the  spherical  wave 
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experiment.  This  occurs  during  the  rebound  phase  of  the  response  and  can  be  seen  in  the  results  shown 
in  Figures  l-10(f)  through  l-lO(h).  A  careful  comparison  of  the  results  shown  in  these  figures  reveals  a 
systematic  difference  between  the  late-time  particle  velocity  histories  recorded  by  gages  in  the  ‘a’  array 
and  those  recorded  by  gages  in  the  ‘b’  array.  A  single  rebound  signal  with  a  duration  of  about  15-17  p^s 
is  observed  on  the  ‘a’  side  of  the  specimen.  In  contrast,  a  two-spike  rebound  signal  is  observed  on  the 
b  side  of  the  specimen.  The  first  spike  has  a  larger  amplitude  and  longer  duration  than  the  second 
spike,  but  when  the  two  spikes  on  the  b  side  of  the  specimen  are  combined,  they  have  the  same 
duration  as  the  single  spike  on  the  ‘a’  side.  We  believe  that  these  systematic  differences  are  closely 
related  to  the  effects  of  cracks  on  wave  propagation.  Specifically,  Ae  lack  of  symmetry  in  the  cracking 
patterns  caused  the  stress  wave  to  become  asymmetric,  and  thus  gages  located  in  different  regions  of  the 
specimen  recorded  different  motions. 

Code  simulations  were  performed  to  aid  in  interpretation  of  the  results  of  the  spherical  wave  experiment. 
The  simulations  were  performed  in  DYNA2D  (Whirley  and  Engelmann,  1992)  using  a  newly  developed, 
computationally  oriented  constitutive  model  based  on  the  evolution  of  damage  by  crack  propagation 
(Simons  et  al.,  1995). 

After  choosing  reasonable  values  for  the  material  parameters  in  the  model  on  the  basis  of  published  data 
for  marble  and  similar  geologic  materials,  we  attempted  to  simulate  the  spherical  wave  experiment. 
However,  with  the  cracking  model  alone,  we  were  unable  to  obtain  a  good  qualitative  agreement 
between  the  measured  and  calculated  particle  velocity  histories.  To  obtain  such  agreement,  we  found  it 
necessary  to  use  a  porous  compaction  model  to  represent  the  material  behavior  in  the  vicinity  of  the 
charge  cavity  (within  a  diameter  of  3.5  cm). 

The  simulated  crack  patterns  in  the  spherical  wave  experiment  are  shown  in  Figure  1-12  along  with  the 
experimentally  observed  crack  patterns.  The  model  reproduced  the  radial  and  circumferential  crack 
networks  observed  in  the  experiment.  The  crack  observed  experimentally  along  the  path  of  the  in  situ 
joint  in  the  rock  was  not  reproduced  in  the  simulation  because  the  joint  itself  was  not  modeled  (i.e., 
initially  isotropic  behavior  was  assumed  in  the  calculations). 

The  model  was  not  quantitatively  calibrated  to  reproduce  all  the  details  of  the  experimental  data.  The 
objective  of  the  simulations  was  to  try  to  provide  a  better  understanding  of  the  mechanisms  that 
contributed  to  the  observed  response.  To  this  end,  our  results  indicate  that  close  to  the  charge  cavity,  the 
behavior  of  the  rock  is  dominated  by  porous  compaction,  whereas  away  from  the  charge,  the  behavior  is 
dominated  by  crack  propagation. 


21 


Circumferential 

Cracks 


1.3  REPORT  ORGANIZATION. 

The  literature  review  is  presented  in  Section  2,  the  uniaxial  strain  experiments  in  Section  3,  and  the 
spherical  wave  experiment  in  Section  4.  The  report  also  includes  three  self-contained  appendices.  The 
first,  Appendix  A,  describes  a  simple  analysis  for  calculating  the  rise  time  of  in-material  gages;  the 
second.  Appendix  B,  summarizes  the  experimental  data  from  all  the  uniaxial  strain  experiments 
performed  during  the  study;  and  the  third.  Appendix  C,  describes  an  analytical  solution  for  the  flow  field 
produced  by  the  reflection  of  compression  waves  at  impedance  mismatched  boundaries.  This  solution 
was  derived  by  M.  Cowperthwaite  (1995)  of  Poulter  Laboratory. 
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SECTION  2 


LITERATURE  REVIEW  OF  THE  SHOCK  RESPONSE  OF  CALCITE 

Investigations  of  the  high  pressure  behavior  of  geologic  materials  under  static  loading  conditions  can  be 
traced  back  to  the  1920s  (Adams  and  Williamson,  1923).  These  early  efforts  focused  primarily  on 
measuring  the  compressibility  of  minerals  of  various  compositions.  Not  until  the  pioneering  work  of 
Bridgman  (1939)  were  phase  transformations  in  geologic  materials  investigated.  The  use  of  shock  wave 
experiments  to  study  the  high  pressure  behavior  of  geologic  media  dates  back  to  the  1960s.  The  early 
work  was  motivated  in  part  by  interest  in  investigating  reactions  and  processes  in  the  earth’s  mantle  that 
occur  at  a  combination  of  pressures  and  temperatures  unattainable  in  static  experiments  (Ringwood, 
1962;  Ahrens  and  Gregson,  1964).  Since  calcite  rocks  such  as  marble  and  limestone  are  among  the  most 
abundant  carbonate  mineralogies  in  the  crust  of  the  earth,  many  investigations  have  been  performed  to 
characterize  their  behavior  under  conditions  of  high  pressure  and  temperature. 

This  section  presents  the  findings  of  a  critical  review  of  existing  studies  concerned  with  the  behavior  of 
calcite  and  calcite-based  geologic  materials  under  high  pressure  and  temperature.  The  material  behavior 
parameters  emphasized  during  this  review  were 

•  Phase  transitions 

•  Loading  Hugoniot 

•  Unloading  adiabat 

•  Effect  of  joints  on  wave  propagation. 

Although  the  review  was  not  exhaustive,  it  allows  us  to  present  here  a  fairly  accurate  account  of  the  state 
of  the  art  of  existing  knowledge  regarding  the  behavior  of  calcite  and  calcite-based  minerals,  especially 
with  regard  to  the  material  behavior  parameters  enumerated  above. 

2.1  PHASE  TRANSITIONS. 

Several  calcite  polymorphs  have  been  identified  in  the  literature.  Calcite  I  is  the  stable  phase  under 
conditions  of  ambient  pressure  and  temperature.  It  has  a  rhombohedral  structure  and  a  crystal  density  of 
2.71  g/cm^.  Aragonite  becomes  the  stable  phase  under  high  pressure.  It  has  an  orthorhombic  structure 
and  a  crystal  density  of  2.93  g/cm^  (Kerley,  1989).  The  phase  diagram  of  calcium  carbonate  rocks  is 
further  complicated  by  the  appearance  of  the  calcite  n  and  calcite  HI  intermediate,  high  pressure  phases 
in  the  pressure  range  from  10  to  25  kbar  (Bridgman,  1939, 1948).  At  even  higher  pressures,  calcite  is 
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known  to  undergo  yet  another  phase  transition  (Adadurov  et  al.,  1960;  Muni  et  al.,  1975;  Grady,  1986); 
this  high  pressure  phase,  known  as  calcite  IV,  occurs  at  pressures  in  excess  of  100  kbar  and  is  usually 
observed  in  impact  experiments  in  the  form  of  a  somewhat  anomalous  rarefaction  shock.  None  of  the 
intermediate  phases  of  calcite  has  a  field  of  stability  under  equilibrium  conditions. 

Kerley  (1989)  developed  a  theoretical  model  for  the  equation  of  state  of  CaCOs  that  includes  phase 
transitions  and  melting.  For  the  solid  phases,  the  model  combines  the  contributions  from  the  cold 
compression  curve  and  lattice  vibrations,  and  for  the  liquid  phase,  the  model  is  constructed  using  fluid 
perturbation  theory.  The  phase  diagram  calculated  using  this  model  is  shown  in  Figure  2-1.  As  noted, 
calcite  IV  is  assumed  to  have  the  same  properties  as  aragonite  except  for  a  shift  in  the  energy. 

The  calcite  H  and  calcite  HI  polymorphs  of  CaCOa  were  first  observed  experimentally  by  Bridgman 
(1939)  at  pressure  levels  of  14.7  and  17.7  kbar,  respectively.  Subsequent  static  measurements  performed 
by  other  investigators  yielded  similar  results.  The  results  of  these  static  studies  are  summarized  in 
Table  2-1.  As  indicated  in  the  table,  the  transition  pressures  and  associated  changes  in  specific  volume 
measured  by  the  various  investigators  are  in  reasonable  agreement  A  noteworthy  exception  is  the  study 
performed  by  Adadurov  et  al.  (1960),  who  identified  three  transitions,  instead  of  two,  in  the  stress  range 
between  0  and  30  kbar.  The  results  of  this  study  are  shown  in  Figure  2-2,  along  with  similar  data 
obtained  by  Singh  and  Kennedy  (1973).  As  shown,  the  second  and  third  transitions  in  the  Adadurov  et 
al.  (1960)  data  agree  favorably  with  the  calcite  I  -  calcite  n  and  calcite  11  -  calcite  III  transitions 
measured  by  Singh  and  Kennedy  (1973).  It  is  not  likely  that  the  first  transition  identified  by  Adadurov 
et  al.  (1960)  is  a  new  transition,  different  from  the  two  identified  by  other  investigators.  It  is  possible 
that  the  first  discontinuity  observed  by  Adadurov  et  al.  (1960)  signifies  the  initiation  of  the  calcite  I  - 
calcite  II  transition  in  small  portions  of  the  specimen  where  local  stress  concentrations  lead  to  stresses 
that  are  higher  than  the  remote  continuum  stresses.  The  second  discontinuity  would  then  signify  the 
completion  of  the  phase  transformation. 

In  some  respects,  this  explanation  is  consistent  with  the  results  of  the  ultrasonic  wave  measurements 
performed  by  Wang  (1966).  Wang  s  results  are  shown  in  Figure  2-3.  The  initiation  of  the  phase 
transition  (at  a  microscopic  level)  is  signified  in  this  figure  by  the  drop  in  compressional  wave  velocity 
as  a  function  of  increasing  pressure.  The  completion  of  the  transition  corresponds  to  the  pressure  level 
where  the  wave  velocity  attains  its  minimum  value.  Figure  2-2  indicates  that,  for  all  the  calcite  rocks 
investigated,  completion  of  the  transition  occurs  within  a  narrow  pressure  range  (15.5  ±  0.5  kbar). 
However,  the  transformation  initiates  at  different  pressures  and  proceeds  at  different  rates  depending  on 
the  morphology  of  the  rock.  In  coarse-grained  rock  (e.g.,  Solenhofen  limestone),  the  transformation 
begins  sooner  and  proceeds  slower  than  in  fine-grained  rock  (e.g.,  Danby  marble).  In  a  single  crystal, 
the  transition  begins  and  ends  abruptly  within  a  narrow  band  of  0.5  kbar.  This  trend  can  be  explained  in 
terms  of  local  stress  fluctuations.  With  increasing  grain  size,  the  material  behaves  less  and  less  like  a 
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TEMPERATURE  (K) 


Figure  2-1 .  Phase  diagram  for  CaC03  including  metastable  phases  II,  III,  and  IV. 

The  calcite  IV  properties  are  the  same  as  for  aragonite,  except  for 
a  shift  in  the  energy  (after  Kerley,  1 989). 
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Table  2-1.  Static  experimental  measurements  of  the  pressures  and  volume 
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Figure  2-3.  Compressional  wave  velocity  as  a  function  of  pressure  for  several  calcite 
rocks.  The  sharp  drop  in  velocity  at  a  pressure  of  about  15  kbar  is 
associated  with  the  calcite  I  -  calcite  II  transition  (Wang,  1966). 
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continuum;  local  stress  gradients  develop  that  lead  to  localized  regions  of  high  stress  concentrations. 
These  regions  of  high  stress  concentrations  transform  before  the  rest  of  the  specimen  and  lead  to  the 
gradual  decrease  in  velocity  observed  in  Figure  2-3. 

Measurements  such  as  those  shown  in  Figure  2-3  may  be  significantly  influenced  by  anisotropy.  This  is 
illustrated  in  Figure  2-4,  where  the  wave  velocity  is  displayed  as  a  function  of  pressure  for  three 
mutually  perpendicular  directions  in  Solenhofen  limestone.  As  shown,  the  pressure  that  corresponds  to 
the  initial  drop  in  velocity  varies  from  one  orientation  to  the  next  by  as  much  as  5  kbar.  The  pressure 
that  corresponds  to  minimum  velocity  varies  between  14  and  15  kbar.  Such  variations  contribute  to 
some  of  the  scatter  in  the  results  presented  in  Table  2-1. 

The  results  shown  in  Figures  2-3  and  2-4  correspond  to  pressures  between  1  and  20  kbar.  In  this 
pressure  range,  Wang  (1966)  observed  only  the  calcite  I  -  calcite  n  transition.  Wang  and  Meltzer  (1973) 
extended  the  pressure  range  by  a  few  kilobars  and  were  able  to  detect  the  calcite  n  -  calcite  HI  transition 
at  a  pressure  of  about  23  kbar,  as  shown  in  Figure  2-5.  The  first  transition  shown  in  this  figure  is  in 
good  agreement  with  the  earlier  results  of  Wang  (1966),  shown  in  Figure  2-3. 

It  has  already  been  shown  that  the  phase  boundaries  among  the  various  polymorphs  of  calcite  arc 
affected  by  temperature  (see  Figure  2-1).  Experimental  measurements  that  verify  this  phenomenon  arc 
shown  in  Figures  2-6  and  2-7.  Figure  2-6  illustrates  this  behavior  in  terms  of  the  shift  in  the  pressure 
that  corresponds  to  minimum  compressional  wave  velocity.  These  results  are  combined  with 
Bridgman’s  (1939)  data  in  Figure  2-7.  The  trends  exhibited  by  the  two  sets  of  data  are  very  similar. 
Quantitative  differences  are  small  enough  to  be  attributed  to  possible  variations  between  the  rock 
specimens  used  in  the  two  different  investigations.  The  results  shown  in  Figure  2-7  are  in  good 
agreement  with  Kerley’s  analytically  derived  phase  diagram  shown  in  Figure  2-1. 

Our  discussion  has  so  far  been  limited  to  static  data,  but  many  wave  propagation  studies  have  been 
undertaken  to  investigate  the  calcite  phase  transformations  under  dynamic  conditions.  Several  of  these 
investigations  are  summarized  in  Table  2-2.  The  dynamic  data  in  this  table  differ  from  their  static 
counterparts  in  that  the  combination  of  high  pressure  and  high  temperature  environments  achievable 
dynamically  are  not  easily  attainable  using  conventional  static  testing  techniques.  For  this  reason,  the 
calcite  IB  -  calcite  IV  transition  observed  in  many  wave  propagation  experiments  was  not  detected  in  the 
static  tests.  This  transition  occurs  in  the  pressure  range  between  95  and  170  kbar.  The  transition 
pressure  of  95  kbar  reported  by  Ahrens  and  Gregson  (1964)  has  not  been  independently  verified  by  other 
investigators  and  appears  to  be  too  low  compared  with  other  measurements. 

Another  difference  between  the  static  and  dynamic  data  is  the  generally  higher  calcite  n  -  calcite  HI 
transition  pressures  recorded  in  impact  experiments.  This  difference  may  be  due  to  the  higher 
temperatures  usually  attained  in  dynamic  experiments.  According  to  the  pressure-temperature  relations 
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Figure  2-4.  Wave  velocity  as  a  function  of  pressure  for  Solenhofen  limestone  in  three 
mutually  perpendicular  directions  (Wang,  1966). 
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Figure  2-5.  Compressional  and  shear  wave  velocities  as  a  function  of  pressure  for 
Oak  Hall  quarry  limestone  (Wang  and  Meltzer,  1973). 
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VELOCITY  (km/s) 


Figure  2-6.  Compressional  wave  velocity  as  a  function  of  pressure  for  Danby  marble 
depicting  the  dependence  of  calcite  I  -  calcite  II  phase  transition  pressure 
on  temperature  (Wang,  1968). 
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TEMPERATURE  (°C) 
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Figure  2-7.  Calcite  I  -  calcite  II  and  calcite  II  -  calcite  III  phase  boundaries  in  pressure- 
temperature  space  for  calcite  rock. 
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Table  2-2.  Dynamic  experimental  measurements  of  the  pressures  associated 
with  the  various  calcite  phase  transitions. 


depicted  in  Figures  2-1  and  2-7,  the  calcite  II  -  calcite  HI  transition  pressure  increases  with  increasing 
temperature.  Furthermore,  the  phase  transition  in  dynamic  experiments  is  usually  observed  during 
unloading  (i.e.,  reverse  transition),  which  means  that  the  specimen  would  have  reached  its  maximum 
temperature  before  the  transition  pressure  was  measured.  If  the  peak  pressure  is  high  and/or  the 
specimen  is  porous,  temperatures  on  the  order  of  few  hundred  degrees  are  easily  attainable  in  an  impact 
experiment,  which  in  turn  leads  to  the  upward  displacement,  in  pressure,  of  the  phase  boundary.  A 
similar  trend  is  not  expected  for  the  calcite  I  -  calcite  n  transition,  simply  because  the  slope  of  the  phase 
boundary  in  the  pressure-temperature  space  is  essentially  flat  in  the  temperature  range  of  interest. 

2.2  LOADING  HUGONIOT. 

Impact  experiments  at  pressures  ranging  from  a  few  kilobars  to  several  hundred  kilobars  have  been 
performed  to  characterize  the  loading  Hugoniot  of  calcium  carbonate  rock  and  to  determine  the  effects 
of  porosity,  saturation,  freezing,  phase  transformations,  and  anisotropy  on  the  dynamic  material 
behavior.  Several  such  investigations  pertinent  to  the  present  study  are  summarized  in  Table  2-3.  The 
Hugoniot  data  from  these  investigations  are  shown  in  Figure  2-8  in  terms  of  stress  versus  volumetric 
strain.  In  this  figure,  the  hollow  symbols  represent  data  for  calcite  in  the  porous  state,  the  solid  symbols 
represent  data  for  calcite  in  the  nonporous  state  (porosity  less  than  0.5%)  and  the  composite  symbols 
(e.g.,  crossed  hollow  circle)  represent  data  for  saturated  and/or  frozen  calcite.  Scatter  in  the  data  is 
evident.  However,  if  the  data  are  separated  into  classes  based  on  porosity  and  saturation,  a  pattern 
begins  to  emerge.  This  observation  is  illustrated  in  Figures  2-9  through  2-11.  Figure  2-9  shows  the  data 
for  nonporous  calcium  carbonate  rock.  The  very  good  correlation  observed  between  the  data  from 
different  investigations  indicates  that  the  high  pressure  equation  of  state  representing  the  loading 
Hugoniot  of  nonporous  calcium  carbonate  rock  is  predictable  with  good  accuracy. 

Figure  2-10  shows  the  results  of  a  study  that  investigated  the  effect  of  saturation  on  the  behavior  of 
Icelandic  spar,  a  nonporous,  single  crystal  calcite  rock.  As  shown,  the  effect  of  saturation  on  the 
behavior  of  the  rock  is  barely  perceptible.  In  a  similar  fashion,  Furnish  (1994),  who  tested  frozen  and 
saturated  low  porosity  Louisville  and  Jacksonville  limestones,  observed  no  systematic  differences 
between  the  behaviors  of  the  saturated  and  the  frozen  samples.  The  results  of  these  two  studies 
collectively  indicate  that,  when  the  rock  is  nonporous,  saturation  and  freezing  have  little  effect  on  the 
Hugoniot. 

Figure  2-11  shows  the  results  of  several  studies  that  investigated  the  effects  of  saturation  and  freezing  on 
the  behavior  of  porous  calcite  rock.  For  each  study  represented  in  the  figure,  the  data  for  the  dry  rock 
are  shown  with  hollow  symbols,  the  data  for  the  saturated  rock  are  shown  with  solid  symbols,  and  the 
data  for  frozen  samples  are  shown  with  composite  symbols  (crossed  hollow  circle).  As  expected,  a 
marked  difference  can  be  observed  between  the  Hugoniot  of  the  dry  samples  and  that  of  the  saturated 
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*The  Danby  (Vermont)  marble  used  in  the  present  investigation  is  the  same  as  that  used  in  the  iarge  scaie  experiments  of  Keough  et  ai.  (1 993). 


Adadurov  et  al.  (1960),  Icelandic  Spar  -  Static 


(JBq>l)  SSHdlS 
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Figure  2-8.  A  coliection  of  Hugoniot  data  for  calcium  carbonate  (CaCOg)  rock  in  the  stress-volumetric  strain  piane. 


400 


in  the  stress-volumetric  strain  plane. 


STRESS  (kbar) 


VOLUMETRIC  STRAIN 


Figure  2-1 0.  Comparison  of  the  ioading  Hugoniots  of  dry  and  saturated  nonporous 
caicium  carbonate  (CaCOg)  rock. 


40 


400 


(and  saturated-frozen)  samples.  Hydration  reactions  in  saturated  calcite,  leading  to  the  formation  of 
ikaite  (CaC03»6H20)  as  suggested  by  Furnish  (1994),  may  be  a  contributing  factor  to  the  observed 
differences. 

Figure  2-11  also  shows  that  the  Hugoniot  data  for  the  dry  samples  from  the  various  investigations  fall 
essentially  on  the  same  curve  as  do  the  data  for  the  saturated  and  frozen  samples.  This  is  because  the 
data  shown  in  the  figure  correspond  to  rock  samples  of  similar  composition  (nearly  pure  CaCX)3)  and 
porosity  (between  14%  and  18%).  An  interesting  observation  based  on  the  data  shown  in  Figure  2-1 1  is 
that  freezing  had  no  detectable  effect  on  the  loading  Hugoniot  of  porous,  saturated  calcite  rock. 

Another  useful  representation  of  the  data  shown  in  Figures  2-8  through  2-1 1  is  in  the  stress-volume 
plane,  as  shown  in  Figure  2-12.  In  this  figure,  as  in  Figure  2-8,  the  hollow  symbols  represent  data  for 
calcite  in  the  porous  state,  the  solid  symbols  represent  data  for  calcite  in  the  nonporous  state  (porosity 
less  than  0.5%),  and  the  composite  symbols  (e.g.,  crossed  hollow  circle)  represent  data  for  saturated 
and/or  frozen  calcite.  The  figure  shows  that,  at  low  pressures  (less  than  about  50  kbar),  significant 
differences  exist  between  the  behaviors  of  rock  of  different  porosities.  With  increasing  pressure,  the 
data  for  both  porous  and  nonporous  rocks  converge  to  the  same  Hugoniot  curve,  which  represents  the 
high  pressure  solid  equation  of  state  of  calcite  IV  -  the  high  pressure  polymorph  of  calcium  carbonate 
rock. 

As  in  the  static  case,  dynamic  data  are  influenced  by  anisotropy.  Ahrens  and  Gregson  (1964)  performed 
plane  wave  experiments  on  Icelandic  spar  single  crystal  specimens.  The  samples  were  oriented  so  that 
the  plane  of  the  shock  was  perpendicular  to  the  x-axis,  y-axis,  z-axis,  and  c-axis  of  the  crystal.  Their 
results  show  a  distinct  difference  in  behavior  among  the  crystals  of  different  orientations  up  to  a  stress  of 
100  kbar.  The  Hugoniot  elastic  limits  (HEL)  measured  for  the  various  crystal  orientations  ranged  from 
18.5  kbar  (c-axis)  to  23.7  kbar  (y-axis). 

2.3  UNLOADING  ADIABAT. 

Generally  speaking,  the  unloading  behavior  of  calcium  carbonate  rock  has  been  less  investigated  than 
the  loading  behavior.  The  recent  study  performed  by  Furnish  (1994)  (see  Table  2-3)  is  the  most 
extensive  investigation  to  date  of  the  unloading  behavior  of  calcite  rocks.  Another  significant  study 
investigating  the  unloading  behavior  of  calcium  carbonate  rock  was  performed  by  Murri  et  al.  (1975). 
Most  of  the  information  provided  in  this  section  is  derived  from  the  results  of  these  two  studies. 

Murri  et  al.  (1975)  performed  plane  wave  experiments  on  Linden  Hall  limestone  and  on  dry  and 
saturated  Indiana  limestone  at  peak  stresses  ranging  between  100  and  300  kbar.  In  all  cases,  rarefaction 
shocks  associated  with  the  calcite  IV  -  calcite  IE,  calcite  in  -  calcite  H,  and  calcite  H  -  calcite  I 


42 


Adadurov  et  al.  (1960),  Icelandic  Spar  -  Static 
Adadurov  et  al.  (1 960),  Icelandic  Spar  -  Dry 
Adadurov  et  al.  (I960),  Icelandic  Spar  -  Saturated 
Ahrens  and  Gregson  (1964),  Spergen  Limestone 
Ahrens  and  Gregson  (1964),  Yule  Marble 
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transformations  were  observed  in  the  experimental  records  during  unloading.  Not  all  of  these 
transformations  were  observed  during  loading.  The  calcite  I  -  calcite  11  transition,  for  example,  was  only 
observed  during  loading  in  the  experiments  on  nonporous  limestone,  so  in  the  porous  case  the  transition 
may  have  been  obscured  by  processes  associated  with  the  compaction  of  the  pores.  It  is  also  noteworthy 
that  the  calcite  11  -  calcite  HI  transition  was  not  observed  during  loading  in  any  of  the  experiments 
conducted  by  Murri  et  al.  (1975),  which  led  them  to  conclude  that  the  kinetics  of  the  calcite  n  -  calcite 
in  transformation  were  too  slow  for  the  transition  to  be  observed  in  a  shock  wave  experiment.  The 
reverse  transformation,  however,  occurred  fast  enough  for  detection.  Porosity  and  water  content  did  not 
seem  to  greatly  affect  the  stress  at  which  the  rarefaction  shocks  associated  with  the  reverse  phase 
transformations  in  calcite  occured. 

Murri  et  al.  (1975)  observed  little  hysteresis  in  the  unloading  behavior  of  nonporous  limestone.  Porous 
limestone,  on  the  other  hand,  exhibited  some  hysteresis  upon  unloading,  in  part  because  of  shock 
heating.  Water-saturated,  porous  limestone  was  intermediate  between  solid  and  dry  porous  limestone. 
This  behavior  is  illustrated  in  Figures  2-13  and  2-14,  which  show  the  unloading  paths  for  dry  and  water- 
saturated  Indiana  limestone,  respectively. 

Furnish  (1994),  who  also  calculated  the  unloading  paths  for  porous  (Indiana  limestone)  and  nonporous 
(Danby  marble)  calcite  up  to  stresses  of  about  130  kbar,  observed  more  hysteresis  in  the  behavior  of  the 
nonporous  samples  than  that  of  the  porous  samples. 

2.4  EFFECT  OF  JOINTS  ON  WAVE  PROPAGATION. 

The  effect  of  joints  on  wave  propagation  in  geologic  material  has  not  been  extensively  studied 
experimentally.  The  two  known  studies  were  performed  by  Gaffney  (1992)  and  by  Gefken  and  Florence 
(1993). 

Gaffney  (1992)  studied  wave  propagation  in  marble  specimens  containing  1-mm-thick  water-filled  and 
ice-filled  artificial  joints  under  uniaxial  strain  conditions  at  stresses  of  10  and  50  kbar.  The  investigation 
was  aimed  at  determining  the  transition  pressure  at  which  the  effect  of  freezing  vanishes  (i.e.,  the 
pressure  beyond  which  water-filled  and  ice-filled  joints  have  the  same  effect  on  wave  propagation). 
Results  in  the  form  of  stress  history  measurements  behind  the  jointed  region  of  the  specimen  indicated 
that  at  10  kbar  ice-filled  and  water-filled  joints  attenuated  the  stress  pulse  by  78%  and  60%,  respectively, 
and  at  50  kbar  ice-filled  and  water-filled  joints  attenuated  the  stress  pulse  by  40%  and  24%,  respectively. 
These  results  indicate  that  the  transition  pressure  under  the  conditions  investigated  in  this  study  is  higher 
than  50  kbar. 

The  second  investigation  of  the  effect  of  joints  on  wave  propagation  in  geological  materials  is  that  of 
Gefken  and  Florence  (1993),  who  used  the  spherical  wave  experimental  technique  to  measure 
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(a)  Particle  velocity  plane  (b)  Stress-specific  volume  plane 

Figure  2-1 3.  Hugoniot  and  release  data  for  dry  Indiana  limestone. 
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(a)  Stress-particle  velocity  plane  (b)  Stress-specific  volume  plane 


Figure  2-14.  Hugoniot  and  release  data  for  water-saturated  Indiana  limestone. 

Heavy  line  is  Hugoniot  and  light  lines  are  release  paths  determined  from  experiment. 
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attenuation  characteristics  of  particle  motion  across  smooth  (polished)  and  rough  (sand-filled)  artificial 

joints  in  Indiana  limestone.  The  study  included  axial  and  radial  velocity  measurements  in  both  dry  and 
saturated  samples. 

Polished  joints  in  dry  limestone  showed  no  indication  of  radial  slip  during  the  passage  of  the  initial 
compression  wave.  Late  time  slip  was  observed  across  the  joint.  Polished  joints  in  saturated  limestone 
showed  evidence  of  instantaneous  slip  as  well  as  late  time  separation  of  the  two  joint  surfaces.  A 
polished  joint  in  dry  limestone  produced  more  transmitted  radial  motion  than  a  similar  joint  in  saturated 
limestone  or  a  rough,  sand-filled  joint  in  either  dry  or  saturated  limestone. 

2£  SYNOPSIS. 

Many  studies  have  been  conducted  over  the  past  30  years  to  characterize  the  behavior  of  calcium 
carbonate  rock  under  dynamic  loadings.  Many  of  the  early  dynamic  wave  propagation  studies 
emphasized  the  characteristics  of  the  loading  Hugoniot  and  to  a  lesser  extent,  the  effects  of  strength, 
porosity,  and  phase  transformations  on  the  material  behavior. 

Static  high  pressure  studies  played  a  significant  role  in  furthering  our  understanding  of  the  behavior  of 
calcium  carbonate  rock.  The  pressure  and  volume  changes  associated  with  the  calcite  I  -  calcite  II  and 
calcite  11  -  calcite  III  phase  transitions  have  been  measured  accurately  in  static  experiments.  These 
transitions  have  since  been  observed  in  shock  wave  propagation  experiments  and  are  now  known  to 
affect  the  shape  of  both  the  loading  Hugoniot  and  the  unloading  adiabat  of  calcite  rock  in  a  predictable 
fashion.  The  loading  and  unloading  behaviors  of  calcite  rock  are  also  affected  by  another  phase 
transformation:  the  calcite  HI  -  calcite  IV  transition,  which  occurs  at  stresses  in  excess  of  100  kbar  and 
has  only  been  observed  in  dynamic  wave  propagation  experiments. 

The  Hugoniot  of  calcium  carbonate  rock  has  been  investigated  for  stresses  up  to  about  500  kbar.  For 
nonporous  rock,  the  Hugoniot  is  well  characterized  and  results  from  numerous  investigations  are  very 
reproducible.  The  Hugoniot  elastic  limit  (HEL)  of  calcite  occurs  in  the  stress  range  from  10  to  20  kbar 
depending,  on  the  porosity  and  graininess  of  the  sample.  The  HEL  is  associated  with  the  first  two 
polymorphic  calcite  transitions  (I  -  n  and  II  -  III). 

The  unloading  adiabat  is  not  too  far  removed  from  the  loading  Hugoniot,  partly  because  of  the  fully 
reversible  nature  of  the  calcite  phase  transitions.  Little  hysteresis  has  been  observed  upon  complete 
unloading  from  a  shock-loaded  state. 

The  effect  of  joints  on  material  response  is  the  least  investigated  aspect  of  the  behavior  of  calcite  rock. 
Our  literature  review  unveiled  only  two  studies  aimed  at  investigating  the  effect  of  joints  on  wave 
propagation.  These  studies  showed  drastic  differences  between  the  response  of  intact  and  jointed 
samples.  This  information,  combined  with  DNA's  interest  in  obtaining  statistical  variations  in  material 
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properties  for  jointed  limestone  and  developing  techniques  for  including  the  observed  effects  in  material 
models  for  ground  shock,  led  to  the  consensus  that  the  experiments  to  be  conducted  during  the  course  of 
the  present  study  should  emphasize  the  effect  of  joints  on  material  behavior. 
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SECTION  3 


UNIAXIAL  STRAIN  EXPERIMENTS 

3.1  MOTIVATIONS  AND  OBJECTIVES. 

Fracture-free  rock  formations  are  rare  in  nature.  For  this  reason,  constitutive  models  for  predicting 
ground  motion  in  the  vicinity  of  a  large  explosion,  such  as  are  required  for  the  HYDRO  PLUS 
methodology,  must  include  the  effect  of  joints  and  inhomogeneities  on  the  material  response. 

The  experiments  described  in  this  section  were  motivated  by  a  lack  of  experimental  data  on  wave 
propagation  in  jointed  rock.  The  experiments  aimed  to  provide  results  that  will  further  our 
understanding  of  the  effect  of  joints  on  wave  propagation  and  provide  a  consistent  set  of  results  that  can 
be  used  to  develop,  calibrate,  or  validate  constitutive  models  suitable  for  use  in  hydrocode  calculations. 

The  study  consisted  of  experimental  and  analytical  investigations.  Experimentally,  uniaxial  strain 
experiments  were  conducted  on  both  intact  and  jointed  calcite  rock  samples.  In  the  experiments,  in¬ 
contact  explosives,  initiated  using  15.24-cm-diameter  (6.0  inch)  plane  wave  lenses,  were  used  to  provide 
shock  loadings  in  the  rock  samples.  Several  specimen  designs  incorporating  different  joint 
configurations  were  investigated.  Experimental  diagnostics  consisted  of  in-material  Lagrangian  particle 
velocity  gages  embedded  between  rock  layers  at  several  serial  locations  within  the  specimen. 
Analytically,  Lagrangian  finite  difference  hydrocode  simulations  were  carried  out  to  help  in  the 
interpretation  of  experimental  results  and  provide  added  insight  into  the  behavior  of  jointed  calcium 
carbonate  rock. 

3.2  EXPERIMENTAL  TECHNIQUES. 

This  section  describes  the  experimental  techniques  used  to  carry  out  the  uniaxial  strain  experiments  on 
intact  and  jointed  calcite  rock  samples.  The  two  geologic  materials  used  in  the  experiments  are 
described  first,  then  the  procedures  used  to  assemble  and  instrument  the  specimens,  and  finally  details 
of  the  experiments  themselves. 

3.2.1  Materials. 

Two  materials,  UTTR  (Utah  Test  Training  Range)  limestone  and  Danby  (Vermont)  marble,  were 
selected  for  use  in  this  investigation.  The  UTTR  limestone  was  selected  because  it  contains  in  situ 
calcite-filled  fractures  that  are  representative  of  naturally  occurring  joints.  Danby  marble  was  selected 
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because  it  is  a  well  characterized  calcite  rock  with  predictable  behavior  in  the  stress  range  of  interest,  so 
that  when  artificial  joints  are  introduced  into  a  specimen,  it  will  be  possible  to  isolate  the  effect  of  joints 
on  the  observed  behavior. 

3.2.1.1.  Danbv  Marble.  The  Danby  marble  (trade  name  “Montclaire”  marble)  used  in  the  present  study 
was  obtained  from  the  marble  quarry  in  Danby,  Vermont.  This  crack-  and  foliation-free  marble  was 
used  in  several  other  investigations  connected  with  the  DISTANT  MOUNTAIN  test  series  (Furnish, 
1994;  Gaffney  and  Smith,  1994;  Keough  et  al.,  1993).  This  Danby  marble  has  a  very  low  porosity  and  is 
composed  of  nearly  pure  calcite.  It  has  a  density  of  2.71  ±  0.9%  g/cm^,  which  is  essentially  the  same  as 
the  2.71  g/cm^  ideal  crystal  density  of  pure  calcite. 

Slabs  of  Danby  marble  were  rough-cut  from  a  stock  with  cross  sectional  dimensions  of  30.5  cm  (12  in.) 
by  20.3  cm  (8  in.).  The  slabs  were  then  precision-ground  fiat  to  within  a  0.05-mm  (0.002-in.)  tolerance 
over  the  entire  surface  of  the  slab. 

3.2.1.2.  UTTR  Limestone.  As  its  name  indicates,  the  UTTR  limestone  is  native  to  the  Utah  Test  and 
Training  Range.  It  was  obtained  from  vertical  drill  holes  near  Candy  Mountain  at  UTTR.  It  is 
composed  primarily  of  calcite,  quartz,  and  dolomite  and  contains  traces  of  other  minerals,  including 
K-feldspar,  pyrite  and  apatite  (Martin  and  Felice,  1992).  The  exact  mineralogical  composition  varies 
with  depth.  The  material  used  in  our  investigation  was  obtained  from  depth  intervals  ranging  between 
24  and  38  m.  From  the  depth  data  and  the  results  of  the  Martin  and  Felice  (1992)  study,  the  following 
mineralogical  and  physical  properties  were  obtained.  The  UTTR  specimen  used  in  HPEOS  (high 
pressure  equation  of  state)  Experiment  1  was  composed  of  62%  calcite,  31%  quartz,  2%  dolomite,  2% 
feldspar,  2%  illite,  and  1%  pyrite.  The  remainder  of  the  UTTR  specimens  used  in  our  investigation  were 
composed  of  approximately  66%  calcite,  28%  quartz,  4%  dolomite,  1%  feldspar,  and  1%  illite.  The 
material  has  a  density  of  2.70  g/cm^,  a  grain  density  of  2.71  g/cm^.  and  a  porosity  of  0.4%.  The  UTTR 
rock  contains  in  situ  calcite-filled  fractxires/joints  ranging  in  thickness  from  less  than  0.02  mm  to  a  few 
millimeters.  The  samples  used  in  the  present  study  contained  a  random  distribution  of  these  in  situ 
joints.  The  joints  in  our  samples  were  less  than  1  mm  thick. 

The  UTTR  limestone  was  supplied  to  SRI  in  the  form  of  15.24-cm-diameter  (6.0-in.)  cores  obtained 
from  a  bore  hole  at  the  Utah  Test  and  Training  Range.  Slabs  were  cut  from  the  limestone  cores  at  an 
angle  of  42°  with  respect  to  the  axis  of  the  core.  This  procedure  resulted  in  ellipse-shaped  slabs  with 
dimensions  of  30.5  and  20.3  cm  along  the  major  and  minor  axes,  respectively.  The  20.3-cm  dimension 
was  the  same  as  the  diameter  of  the  explosives.  The  30.5-cm  dimension  was  larger  than  the  diameter  of 
the  explosives  and  provided  the  shielding  necessary  to  protect  the  particle  velocity  gages  from  the  blast 
generated  by  the  explosives.  The  slabs  were  precision-ground  fiat  to  within  a  0.05-mm  (0.002-in.) 
tolerance  over  the  entire  surface  of  the  slab. 
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3.2.2  Specimen  Conflgurations  and  Assembly. 

The  specimen  configurations  investigated  during  this  study  are  shown  in  Figure  3-1.  They  included 

•  Intact  Danby  marble  specimens  with  no  joints 

•  Artificially  jointed  Danby  marble  specimens 

•  U'lTR  limestone  specimens  with  in  situ  joints. 

The  same  basic  experimental  design  shown  in  Figure  3-1  was  used  in  all  the  uniaxial  strain  experiments 
except  HPEOS  Experiment  1,  which  is  described  in  a  later  section.  As  shown  in  Figure  3- 1(a),  two 
particle  velocity  gages  on  each  of  four  gage  planes  were  embedded  within  the  specimen.  The  gages  on 
the  first  gage  plane  (closest  to  the  top  surface  of  the  specimen)  were  used  to  record  the  oncoming  stress 
wave  generated  by  the  explosives  and  thereby  provide  a  basis  for  comparing  the  results  of  different 
experiments  with  different  joint  configurations.  In  experiments  on  jointed  rock  samples,  this  gage  was 
also  used  to  record  the  reflected  stress  pulse  from  the  rock-joint  interface.  The  remaining  gages  were 
used  to  record  the  stress  wave  profile  after  the  wave  propagated  through  the  jointed  section  of  the 
specimen.  The  results  from  several  experiments  were  compared  to  determine  the  effect  of  joints  of 
various  configurations  on  the  propagation  of  stress  waves  in  rock. 

Most  experiments  were  conducted  at  a  peak  stress  of  about  120  kbar.  In  each  experiment,  the  shock 
loading  was  provided  using  a  15.24-cm-diameter  (6.0-in.)  plane  wave  lens  and  an  explosive  pad  of  a 
type  and  thickness  appropriate  to  achieve  the  required  loading  conditions  for  the  experiment  (a  5.08-cm- 
thick  BARATOL  pad  was  used  in  most  experiments). 

3.2.2.I.  Artificial  .Toints.  Several  artificial  joint  configurations  were  investigated  during  this  study. 
These  configurations  included  variations  of  surface  texture  (smooth  versus  rough),  thickness  (closed 
versus  sand-filled),  number  (single  versus  triple),  and  orientation  (normal  versus  inclined  with  respect  to 
the  direction  of  shock  propagation).  The  smooth  joint  surfaces  were  precision  ground  as  described 
above.  The  rough  joint  surfaces  were  prepared  by  first  grinding  the  surfaces  flat,  then  inscribing  them 
with  regularly  spaced  grooves  as  shown  in  Figure  3-2.  As  indicated  in  the  figure,  rough  joints  were 
filled  with  a  2-mm-thick  layer  of  Reid-Bedford  sand,  the  properties  of  which  are  described  by  Phillips 
(1986a,  1986b).  Some  sand-filled  smooth  joints  were  also  investigated  during  the  present  study. 

Specimens  with  sand-filled  joints  were  constructed  as  follows.  First,  the  specimen  was  assembled  with 
2.0-mm-thick  shims  around  the  periphery  of  the  joint  surface.  After  the  specimen  construction  was 
complete  and  the  gages  were  properly  wired,  sand  was  pored  into  the  joint(s)  through  egress  holes  along 
the  side  of  the  specimen.  The  sand  was  compacted  by  repeated  tamping  of  the  specimen.  Once  the  sand 
could  no  longer  be  compacted  using  this  method,  the  egress  holes  were  sealed  with  epoxy. 
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Rock  Layer  No.  2a 

No.  2b 

Artificial  Joint 

(b)  Modification  of  the  basic  configuration  to  include  a  single,  normal  joint. 


(c)  Modification  of  the  basic  configuration  to  include  a  triple,  normal  joint. 


Rock  Layer  No.  No.  2b  ^^Artificial  Joint 


(d)  Modification  of  the  basic  configuration  to  include  a  single,  inclined  joint. 


Figure  3-1 .  Configurations  of  the  uniaxial  strain  experiments. 


Transmitted 
Joint  Surface 


Figure  3-2.  Geometry  of  the  surfaces  of  a  rough  joint. 
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3.2.3  Instrumentation. 


Two  types  of  in-material  Lagrangian  gage  systems  were  used  to  record  either  stress  or  particle  velocity 
histories  in  the  impacted  geological  specimens.  Ideally,  a  Lagrangian  gage  monitors  and  records  the 
flow  field  at  some  location  within  the  specimen  without  perturbing  the  motion  of  the  specimen.  In 
practice,  however,  this  is  not  the  case.  Because  of  the  finite  thickness  of  the  gage  package  and  the 
impedance  mismatch  between  the  specimen  material  and  the  gage  plane,  the  flow  field  in  the  vicinity  of 
the  gage  is  disturbed,  and  a  finite  rise  time  is  required  to  establish  a  state  of  equilibrium  at  the  specimen- 
gage  interface.  An  analytical  method  was  devised  during  this  study  to  estimate  the  rise  time  of  in¬ 
material  Lagrangian  gages.  This  method,  which  is  described  in  Appendix  A,  shows  that  the  rise  time  of 
the  gages  used  in  our  investigation  is  less  than  50  ns.  This  is  well  within  the  resolution  for  the 
measurement. 

3.2.3.I.  Stress  Gage  System.  The  stress  gage  system  was  used  only  in  HPEOS  Experiment  1,  which 
was  a  relatively  simple  experiment  designed  primarily  to  investigate  the  effect  of  quartz  content  in  the 
UTTR  rock  on  the  signal-to-noise  ratio  of  transient  measurements  made  using  in-material  stress  or 
particle  velocity  gages.  The  configuration  for  this  experiment  is  shown  in  Figiue  3-3.  As  shown,  the 
specimen  consists  of  a  four-terminal  manganin  gage  sandwiched  between  two  slabs  of  UTTR  limestone. 
The  gage  leads  exit  the  specimen  at  a  45°  angle  so  that  they  are  shielded  from  the  blast  generated  by  the 
explosives  and  a  longer  recording  time  is  ensured.  The  explosives  used  in  this  experiment  consisted  of  a 
10.2-cm-diameter  (4-in.)  plane  wave  lens  and  a  5.08-cm-thick  (2-in.)  pad  of  the  high  explosive  PBX. 

The  four-lead  manganin  gage  grid  shown  in  Figure  3-3  was  photo-etched  from  a  0.05-mm-thick 
manganin  foil.  This  type  of  stress  gage  is  based  on  the  principle  of  piezoresistance.  The  resistance  of 
the  gage,  initially  about  0.1  Q,  changes  predictably  when  a  pressure  disturbance  like  a  stress  wave 
passes  through  the  gage  element.  The  voltage  associated  with  this  change  in  resistance  is  recorded  using 
a  digital  oscilloscope  and  is  then  converted  to  stress  using  well  established  calibration  relationships. 

In  reducing  the  voltage  data,  it  is  assumed  that 


AR  AV 
Ro  ■  Vo 


(3-1) 


where  Ro  is  the  initial  resistance  of  the  gage,  AR  is  the  change  in  resistance,  Vq  is  the  initial  voltage,  and 
AV  is  the  change  in  voltage.  This  assumption  carries  the  implication  that  the  effects  of  changes  in  the 
resistance  of  the  gage  element  and  leads  and  decay  of  the  capacitor  voltage  (power  supply)  are 
negligible.  In  reality,  these  effects  account  for  1%,  or  less,  of  the  recorded  signal. 
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(b)  Top  view  of  the  gage  plane 


Figure  3-3.  Specimen  configuration  and  dimensions  for  HPEOS 
Experiment  1  in  UTTR  limestone. 


55 


In  the  stress  range  of  interest  here,  and  assuming  uniaxial  strain  behavior,  the  relationship  between  the 
change  in  resistance  (or  voltage)  and  the  stress  in  the  manganin  gage  is  given  by  the  equation 


a  (kbar)  =  a 


AR(fl)^ 

Ro(fl); 


(3-2) 


where  a  is  the  stress  and  a  is  an  experimentally  determined  proportionality  constant  with  a  value  of 
350kbar/n/a  (DeCarli,  1976). 

3.2.3  Particle  Velocity  Gage  System.  The  magnetic  particle  velocity  gages  used  throughout  this 
study  are  based  on  the  principle  that  a  conductor  (the  particle  velocity  gage)  moving  in  the  direction 
normal  to  a  uniform  quasi-stadc  magnetic  field  generates  a  voltage  proportional  to  its  length  and 
velocity  and  the  strength  of  the  magnetic  field.  Mathematically,  this  proportionality  is  expressed 
through  the  relation 

u(t)  =  -^xl04  (3-3) 


where  t  is  time,  u(t)  is  particle  velocity  in  millimeters  per  microsecond,  V(t)  is  gage  voltage  in  volts,  t  is 
the  center-to-center  active  gage  length  in  millimeters,  and  B  is  the  magnetic  field  strength  in  gauss. 

The  uniform  magnetic  field  was  provided  using  Helmholtz  coils.  A  schematic  representation  of  a 
specimen  situated  within  a  Helmholtz  coil  assembly  is  shown  in  Figure  3-4.  Each  coil  consisted  of 
40  circular  turns  of  No.  8  copper  wire.  The  coil  radius,  which  was  equal  to  the  separation  distance 
between  the  coils,  was  about  45.72-cm  (18-in.).  The  coil  assembly  was  pulsed  with  a  constant  current  of 
about  350  amperes  a  few  hundred  microseconds  before  the  experiment.  The  current  in  the  coils  was 
measured  using  a  pulse  current  transformer  (Pearson  Electronics,  Inc.,  Palo  Alto,  CA  94303).  The 
magnitude  of  the  uniform  magnetic  field  was  calculated  from  the  relation 

JcNI 

B=  0.286™  (3-4) 


where  B  is  the  field  strength  in  gauss,  N  is  the  number  of  turns  in  one  coil,  I  is  the  current  through  the 
coils  in  amperes,  and  r  is  the  coil  radius  in  centimeters. 

Two  types  of  particle  velocity  gages  were  used  in  this  investigation:  SRI-manufactured  gages 
and  prefabricated  gages  manufactured  by  TAYCO  Engineering  (P.O.  Box  6034,  Cypress,  CA  90630). 
A  detailed  engineering  drawing  of  the  TAYCO  particle  velocity  gages  is  shown  in  Figure  3-5.  As 
shown,  the  gages  were  fabricated  in  pairs,  with  one  gage  made  of  copper  and  the  other  of  Inconel.  Ttie 
gages  were  laid  out  on  a  0.05-mm-thick  Teflon  backing.  A  gage  package  such  as  the  one  shown  in 
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Helmholtz 

Coils 


Figure  3-4.  Schematic  diagram  of  a  typical  uniaxial  strain  experiment 
instrumented  with  particle  velocity  gages. 
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Figure  3-5.  Engineering  drawing  showing  the  specifications  and  dimensions  of  a  prefabricated  particle  velocity  gage. 


Figure  3-5  was  embedded  in  each  gage  plane  to  provide  a  duplicate  recording  of  the  pardcle  velocity  at 
the  gage  plane.  The  active  elements  of  the  prefabricated  gages  each  had  an  area  of  12.7  by  3.18  mm. 
The  SRI  gages  had  the  same  dimensions  as  the  prefabricated  gages  and  had  Inconel  active  elements 
connected  to  copper  leads.  Unlike  the  prefabricated  gages,  the  SRI  gages  were  laid  out  directly  on  the 
specimen  without  any  backing  materials. 

33  EXPERIMENTAL  RESULTS. 

The  measurements  from  the  uniaxial  strain,  in-contact  explosive  experiments  are  presented  in 
Appendix  B.  The  appendix  describes  the  experimental  conditions  and  specimen  geometries  and 
summarizes  the  raw  data  and  reduced  experimental  results.  For  easy  reference.  Appendix  B  includes  a 
tabular  summary  of  all  the  uniaxial  strain  experiments  (Table  B-1). 

In  general,  the  experimental  results  were  self-consistent  and  revealed  important  information  about  phase 
transformations  in  calcium  carbonate  rock  as  well  as  the  effect  of  joints  on  wave  propagation.  The 
results  were  repeatable.  In  each  of  the  uniaxial  strain  experiments,  two  gages  were  used  to  provide 
duplicate  measurements  of  the  particle  velocity  histories  at  each  gage  plane.  In  general,  the  agreement 
between  the  gages  on  a  given  gage  plane  was  excellent,  as  can  be  seen  by  a  close  examination  of  the 
experimental  records  in  Appendix  B. 

Two  types  of  gages  were  used  to  measure  particle  velocity  histories:  SRI-manufactured  gages  and 
prefabricated  gages.  To  ensure  that  both  tj^es  of  gages  yielded  identical  results,  duplicate 
measurements  were  taken  during  two  experiments  (HPEOS  Experiments  6  and  7)  and  the  measured 
histories  were  compared.  The  two  types  of  gages  yielded  identical  results. 

HPEOS  Experiment  1  was  conducted  to  determine  the  effect  of  the  quartz  content  of  the  UTTR 
limestone  on  the  signal-to-noise  ratio  in  the  recorded  data.  The  results  of  this  experiment  (shown  in 
Figures  3-6  and  3-7)  indicate  that  quartz  content  has  little  effect  on  the  recorded  signal.  Figure  3-6 
shows  the  voltage  history  recorded  during  the  experiment.  The  voltage  in  the  gage  element  is  zero 
initialy,  but  increases  to  about  IV  at  t  =  0.5  fis.  This  is  the  exitation  voltage,  and  it  is  provided  using  a 

constant  current  power  supply.  The  voltage  in  the  gage  remains  constant  until  the  arrival  of  the  stress 
wave  at  the  gage  plane  approximatly  20  p.s  after  the  gage  is  pulsed. 

The  voltage  history  shown  in  Figure  3-6  was  converted  to  a  stress  history  by  using  Equation  (3-2),  and 
the  results  are  plotted  in  Figure  3-7.  The  origin  of  the  time  axis  in  this  figure  is  chosen  arbitrarily  and 
has  no  bearing  on  the  results.  The  wave  structure  depicted  in  the  figure  is  relatively  simple:  a 
compressional  loading  wave  that  shocks  the  material  up  to  a  peak  stress  of  260  kbar,  a  gentle  unloading 
wave  that  causes  a  50-kbar  decrease  in  peak  stress,  and  a  main  release  wave  that  unloads  the  rock  to  a 
stress-free  state.  The  initial  shock  wave  originates  in  the  explosive  as  does  the  gentle,  Taylor  unloading 
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Figure  3-6.  Manganin  gage  record  in  HPEOS  Experiment  1 . 
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Figure  3-7.  Stress  history  at  the  gage  plane  in  HPEOS  Experiment  1 . 
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wave.  The  main  release  wave  is  a  backward  facing  wave  that  originates  at  the  back  free  surface  of  the 
specimen.  The  pulse  duration  of  1.8  ^is  is  equivalent  to  the  time  required  for  the  stress  wave  to 
propagate  from  the  gage  plane  to  the  back  surface  of  the  specimen  and  back  again  to  the  gage  plane. 

This  dme  interval  is  less  than  the  time  required  for  the  release  waves  from  the  lateral  specimen 
boundaries  to  reach  the  gage  plane  (approximatly  3  \is),  and  thus  a  uniaxial  strain  state  is  ensured  in  the 

specimen  for  the  duration  of  the  experiment. 

HPEOS  Experiment  2  was  instrumented  as  shown  in  Figure  3-1.  In  addition,  a  particle  velocity  gage 
was  installed  at  the  back  surface  of  the  specimen  to  attempt  to  measure  the  free  surface  velocity  of  the 
sample.  We  believe  that  the  data  recorded  by  this  gage  do  not  reflect  the  real  free  surface  velocity  of  the 
specimen,  for  reasons  illustrated  in  Figure  3-8.  The  figure  shows  a  distance-time  diagram  and  a 
pressure-particle  velocity  diagram  for  the  first  few  wave  reverberations  resulting  from  the  interaction 
between  the  free  surface  of  the  specimen  and  the  gage  package.  Solid  lines  in  Figure  3-8(a)  represent 
compression  waves,  whereas  dashed  lines  represent  rarefaction  waves.  In  reality,  each  of  the  dashed 
lines  represents  the  “toe”  of  a  dispersive  rarefaction  fan,  the  remainder  of  which  is  omitted  from  the 
figure  for  clarity.  Figure  3-8  shows  three  stress  states  in  the  gage  package.  State  “1”  is  a  compressive 
stress  state  that  results  from  the  interaction  of  the  main  pressure  wave  with  the  gage/specimen  interface. 
State  “2”  is  a  stress-free  state  created  by  the  reflection  of  the  stress  wave  from  the  free  surface  of  the 
gage  package.  State  “3”  is  a  tensile  stress  state  due  to  the  second  interaction  between  the  specimen  and 
the  gage  package.  It  is  this  latter  stress  state  that  causes  the  free  surface  velocity  measurement  to  be 
erroneous.  Because  of  the  tensile  stress  generated  at  the  gage/sample  interface,  the  gage  separates  from 
the  back  surface  of  the  specimen  and  moves  at  a  velocity  higher  than  the  velocity  of  the  specimen.  This 
observation  is  supported  by  the  experimental  data  of  HPEOS  Experiment  2,  which  show  the  peak 
velocity  recorded  by  the  free  surface  gage  to  be  more  than  50%  larger  than  the  velocity  recorded  3.3-mm 
away  in  the  interior  of  the  sample.  Ideally,  these  two  gages  should  reach  approximately  the  same 
velocity.  Owing  to  the  failure  at  the  gage/sample  interface,  the  free  surface  gage  was  omitted  from  later 
experiments. 

The  voltage  output  of  the  particle  velocity  gages  in  HPEOS  Experiments  3  and  4  could  not  be  converted 
to  particle  velocity  histories  because  of  failure  of  the  pulse  current  transformer.  This  failure  was 
corrected  and  did  not  affect  subsequent  experiments.  Given  the  consistency  and  good  reproducibility  of 
the  results,  and  the  fact  that  most  of  the  uniaxial  strain  experiments  were  conducted  under  the  same 
boundary  and  loading  conditions,  the  particle  velocity  histories  for  Experiments  3  and  4  can  be 
calculated  with  reasonable  accuracy  (±5%)  by  using  the  calibration  values  of  other  similar  experiments. 
However,  such  calculations  were  not  attempted. 
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(pressure  is  positive  in  compression). 


3.4  LAGRANGIAN  ANALYSIS. 

Since  its  introduction  by  Fowles  and  Williams  (1970),  the  Lagrangian  analysis  method  has  had  a  long 
history  of  improvements  (Cowperthwaite  and  Williams,  1971;  Grady,  1973),  which  culminated  in  the 
development  of  the  GUINSY  computer  code  developed  by  Seaman  (1987)  in  our  laboratory  and 
considered  to  be  the  industry  standard. 

Originally,  we  had  intended  to  use  the  Lagrangian  analysis  code  GUINSY  to  obtain  loading  and 
unloading  paths  for  the  rock  samples  investigated  during  our  study.  However,  this  analysis  could  not  be 
performed  owing  to  inherent  deficiencies  in  the  present  version  of  GUINSY  in  handling  wave 
interactions  resulting  from  reflections  from  impedance-mismatched  boundaries,  as  discussed  below. 

In  the  uniaxial  strain  experiments,  the  compressive  wave  generated  by  the  explosive  travels  through  the 
gaged  region  of  the  specimen,  reaches  the  free  surface,  then  reflects  back  into  the  specimen  as  a  release 
wave  that  unloads  the  specimen  to  a  stress-free  state.  The  forward-traveling  compression  wave  and  the 
backward-traveling  release  wave  interact  near  the  free  surface  of  the  specimen  to  create  a  region  of 
mixed  flow.  Wave  interactions  in  this  region  require  the  use  of  complicated  algorithms  different  from 
the  ones  used  in  routine  GUINSY  applications.  A  method  of  analysis  incorporating  these  algorithms 
was  programmed  into  GUINSY  and  tested  on  our  data.  Because  of  the  complex  nature  of  this  method,  it 
had  been  speculated  when  the  method  was  developed  that  further  developments  might  be  required  before 
it  could  be  used  reliably  (see  Seaman,  1987,  p.  12-14). 

When  an  attempt  was  made  to  use  GUINSY  to  analyze  the  data  from  the  present  experiments, 
deficiencies  were  indeed  discovered  in  the  code’s  ability  to  handle  wave  interactions.  Several  steps  have 
since  been  taken  to  fix  the  problem.  First,  the  old  GUINSY  algorithms  have  been  revised  and  new 
algorithms  have  been  added  to  the  code.  Second,  a  closed-form  analytic  solution  of  a  problem  that 
involves  wave  reflections  from  an  impedance-mismatched  boundary  was  obtained  to  be  used  as  a  test 
case  of  the  new  code  to  ensure  proper  operations  within  the  code.  This  analytic  solution  is  described  in 
detail  by  Cowperthwaite  (1995)  in  the  technical  paper  included  in  Appendix  C. 

We  are  currently  debugging  the  new  algorithms,  and  at  the  time  this  report  was  prepared  the  GUINSY 
code  was  not  robust  enough  to  be  used  reliably  in  the  analysis  of  the  experimental  data.  Instead, 
numerical  simulations  of  the  experiments  were  performed  using  the  Lagrangian  finite  difference 
hydrocode  PUFF  (Seaman  and  Curran,  1978).  These  simulations  are  described  below. 
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3^  CODE  SIMULATIONS. 


Numerical  simulations  were  carried  out  using  the  finite  difference  code  PUFF  to  help  in  the 
interpretation  of  experimental  results  and  provide  added  insight  into  the  behavior  of  jointed  calcium 
carbonate  rock.  This  section  describes  the  constitutive  models  used  in  the  simulations.  It  also  shows  a 
verification  analysis  in  which  PUFF  was  used  to  simulate  an  experiment  similar  to  the  ones  performed 
during  the  present  investigation.  Simulation  results  for  jointed  samples  will  be  described  and  compared 
with  experimental  data  as  appropriate  in  a  later  section. 

In  the  simulations,  the  rock  behavior  was  described  using  a  rate-independent,  elastic-plastic  model  for 
deviatoric  stress  calculations  and  a  Mie-Griineisen  equation  of  state  for  pressure  calculations.  In  PUFF, 
the  Mie-Griineisen  equation  takes  the  following  form: 

P  =  (Cp  +  Dp2  +  s^i3)  ^1 .  rpE  (3-5) 

where  P  is  the  pressure,  F  is  the  Griineisen  coefficient,  p  is  the  density,  E  is  the  internal  energy, 

(Cp,  +  Dp2  +  is  the  Hugoniot  pressure,  C,  D  and  S  are  coefficients  with  moduli  units,  and  |X  is  the 
volumetric  strain  given  by  the  relation 


=  1  (3-6) 

Po 

where  po  is  the  initial  density.  The  model  parameters  were  determined  from  the  collection  of  existing 
data  for  nonporous  calcite  rock  encountered  in  the  literature  review.  These  parameters  are  summarized 
in  Table  3-1,  and  a  comparison  between  experimental  data  for  nonporous  calcite  and  the  Hugoniot 
equation  of  state  used  in  the  code  calculations  is  shown  in  Figure  3-9. 

When  the  response  of  a  specimen  with  a  sand-filled  joint  was  simulated,  the  behavior  of  the  sand  was 
modeled  using  a  porous  equation  of  state.  The  sand  was  assumed  to  have  an  initial  density  of 
1.74  g/cm^,  an  initial  bulk  modulus  of  28  kbar,  and  a  shear  modulus  of  10  kbar  and  to  reach  its  solid 
density  at  a  pressure  of  20  kbar  and  thereafter  behave  like  solid  quartz. 

Before  we  simulated  the  experiments  of  the  present  investigation,  we  carried  out  a  verification  analysis 
in  which  an  experiment  similar  to  the  ones  described  here  was  simulated  and  the  results  compared  to 
experimental  data  [Murri  et  al.  (1975),  Experiment  1883-53].  The  experimental  configuration  is  shown 
in  Figure  3-10.  The  material  used  in  the  experiment  was  Linden  Hall  limestone  (<0.5%  porosity).  The 
in-contact  explosives  used  for  loading  consisted  of  a  20.32-cm-diameter  (8  in.)  plane  wave  lens  and  a 
15.24-cm-thick  (6-in.)  pad  of  Baratol  explosives  [see  Murri  et  al.  (1975)  for  more  details]. 
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Table  3-1.  Parameters  used  to  describe  the  behavior 
of  marble  In  the  PUFF  finite  difference 
simulations. 


Parameter _ Value  Units 


Initial  densit/ 

2.7 

g/cm® 

Yield  stress 

10.0 

kbar 

Shear  modulus 

400.0 

kbar 

Hugoniot  parameters 

C  (bulk  modulus) 

191.7 

kbar 

D 

1732.0 

kbar 

S 

-742.5 

kbar 

Griineisen  coefficient 

1.59 
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STRESS  (kbar) 


Model  Used  in  PUFF  Calculations 
Adadurov  et  al  (1960)  -  Icelandic  Spar 
Gaffney  and  Smith  (1994)  -  Danby  Marble 
Muni  et  al.  (1975)  -  Linden  Hall  Limestone 
Grady  et  al.  (1978)  -  Oak  Hall  Limestone 
Grady  et  al.  (1978)  -  Vermont  Marble 
Grady  (1986)  -  Icelandic  Spar 


□  y  V 


VOLUMETRIC  STRAIN 


Figure  3-9.  Comparison  between  the  equation  of  state  used  in  the  PUFF 
simulations  and  experimental  data  for  nonporous  calcite. 
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Figure  3-10.  Target  configuration  and  dimensions  for  experiment  1 883-53  on 
Linden  Hall  limestone  (Muni  et  al.,  1 975.) 
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The  experimental  and  simulated  particle  velocity  histories  are  compared  in  Figure  3-11.  The  correlation 
between  the  experimental  and  calculated  records  is  rather  close.  The  overall  shape,  first  particle  velocity 
plateau  (corresponding  to  peak  stress),  and  peak  particle  velocities  (corresponding  approximately  to  the 
firee  surface  velocity)  are  all  reproduced  by  the  code  calculations.  The  "kink"  in  each  of  the 
experimental  velocity  histories  just  prior  to  the  peak  velocity  plateau  is  associated  with  a  reverse  phase 
transformation.  This  feature  was  not  reproduced  in  the  code  calculations  because  phase  transitions  were 
not  incorporated  in  the  material  description  in  the  hydrocode. 

The  close  correlation  between  the  simulated  velocity  histories  and  the  experimental  data  in  the  example 
tested  lends  credence  to  the  simulation  results  and  justifies  the  use  of  PUFF  to  help  in  the  interpretation 
of  experimental  results  and  to  provide  added  insight  into  the  effect  of  joints  on  wave  propagation. 

3.6  DISCUSSION. 

3.6.1  Hugoniot  Elastic  Limit  and  Phase  Transitions. 

A  precursor  to  the  main  compression  wave  was  observed  in  the  uniaxial  strain  experiments  on  both 
Danby  marble  and  UTTR  limestone.  The  magnitude  of  the  precursor  was  determined  for  HPEOS 
Experiments  2, 6, 7, 8,  and  16  by  using  the  conservation  of  momenmm  equation  across  the  shock  front, 
namely. 


-  Oo  =  Po  (U-  Uo)(u  -  Uo)  (3-7) 

where  o  and  u  are  the  stress  and  particle  velocity  behind  the  shock  front;  po,  Oq,  and  Uo  are  the  density, 
stress,  and  particle  velocity  ahead  of  the  shock  front;  and  U  is  the  shock  velocity.  Ahead  of  the 
precmsor,  both  Oq  and  Uq  are  zero,  and  Equation  (3-7)  reduces  to 


c  =  poUu.  (3-8) 

Equation  (3-8)  was  used  to  determine  the  precursor  stress  for  Danby  marble  and  UTTR  limestone  by 
Using  the  experimentally  recorded  values  of  particle  velocity  and  shock  velocity.  The  particle  velocity 
along  the  precursor  was  obtained  directly  from  the  experimental  records.  The  precursor  wave  velocity 
was  determined  based  on  the  time  of  arrival  of  the  elastic  wave  at  two  subsequent  gage  locations  and  the 
separation  distance  between  the  gage  planes. 

The  average  precursor  stress  for  Danby  marble  was  17.1  kbar,  based  on  calculated  values  of  16.8, 16.3, 
19.0,  and  16.2  kbar  for  HPEOS  Experiments  2, 6, 7  and  8,  respectively.  These  results  are  in  good 
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PARTICLE  VELOCITY  (mm/|xs)  PARTICLE  VELOCITY  (mm/M,s) 


TIME  (ns) 

(b)  PUFF  simulation 


Figure  3-1 1 .  Comparison  of  experimental  and  calculated 
particle  velocity  histories  for  Murri  et  al. 

(1 975)  Experiment  1 883-53. 
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agreement  with  the  Danby  marble  results  of  Keough  et  al.  (1993),  who  reported  precursor  stress  values 
between  14  and  16  kbar,  Furnish  (1994),  who  reported  precursor  stress  values  between  13  and  15  kbar, 
and  Gaffney  et  al.  (1994),  who  reported  precursor  stress  values  between  14  and  16  kbar. 

The  precursor  stress  for  UTTR  limestone  was  14.3  kbar,  which  was  calculated  based  on  the  results  of 
HPEOS  Experiment  16.  This  stress  is  slightly  higher  than  the  1 1-14  kbar  precursor  stresses  reported  by 
Furnish  (1994)  for  UTTR  limestone. 

Figure  3-12  shows  an  average  of  the  two  particle  velocity  histories  recorded  at  the  second  gage  plane  of 
HPEOS  Experiment  2.  The  prominent  features  of  the  experimental  data  for  intact  marble  are  indicated 
in  the  figure  by  the  letters  A,  A',  B  and  B'.  The  precursor  observed  in  our  experiments  and  indicated  by 
Point  A  in  Figure  3-12  is  most  likely  the  result  of  processes  associated  with  the  calcite  I  -  calcite  n  and 
calcite  II  -  calcite  IB  transitions,  not  plastic  yielding.  This  observation  is  supported  by  our  unloading 
data.  If  the  precursor  were  due  to  yielding,  reverse  yielding  would  be  expected  to  occur  early  during 
unloading  (for  elastic/perfectly  plastic  materials,  reverse  yielding  occurs  when  the  stress  reaches  a  value 
equal  to  the  peak  stress  minus  twice  the  yield  stress).  However,  this  is  not  observed  in  our  results. 
Instead,  the  signature  of  a  reverse  phase  transition  is  observed  in  the  form  of  a  rarefaction  shock  that 
occurs  late  during  unloading,  at  a  stress  level  approximately  equal  to  the  precursor  stress.  This  feature 
(Point  A'  in  Figure  3-10)  indicates  that  the  precursor  observed  in  our  experiments  is  associated  with  the 
calcite  phase  transformations. 

Another  interesting  feature  of  the  particle  velocity  history  shown  in  Figure  3-12  is  an  apparent  high- 
pressure  phase  transition  at  Point  B  and  a  corresponding  reverse  phase  ttansition  at  Point  B’.  Assuming 
a  two  wave  structure,  and  making  use  of  Equation  (3-7)  and  the  mass  conservation  equation  across  a 
shock  front. 


P__i  u-Uq 
Po  U-Uo 


(3-9) 


the  pressure  at  Point  B  can  be  calculated.  For  HPEOS  Experiment  2,  this  pressure  is  1 18.3  kbar.  High 
pressure  phase  transitions  in  calcium  carbonate  rocks  have  been  observed  in  this  pressure  range  by  Murri 
et  al.  (1975),  Ahrens  and  Giegson  (1964),  and  Adadurov  et  al.  (1960)  (see  Table  2-2),  findings  providing 
further  evidence  that  the  wave  structure  at  points  B  and  B'  is  associated  with  this  high  pressure  phase 
transition  in  calcite  rock. 


3.6.2  Effect  of  Joints  on  Wave  Propagation. 

Experiments  were  performed  to  determine  the  effects  of  in  situ  and  artificial  joints  on  wave  propagation 
in  CaC03  rock.  The  in  situ  joints  in  UTTR  limestone  did  not  have  a  measurable  effect,  as  evidenced  by 
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the  similarity  between  the  particle  velocity  histories  recorded  in  the  UTTR  limestone  experiment 
(HPEOS  Experiment  16)  and  their  counterparts  in  the  intact  marble  experiments  (HPEOS  Experiments 
2, 6,  and  7).  On  the  other  hand,  artificial  joints  in  marble  had  a  significant  effect  on  the  observed 
response.  To  better  understand  the  effect  of  artificial  joints  on  wave  propagation,  we  performed  a  series 
of  hydrocode  simulations  using  the  finite  difference  wave  propagation  hydiocode  SRI  PUFF  (Seaman 
and  Curran,  1978).  First,  we  simulated  the  response  of  an  intact  marble  specimen  having  the 
configuration  shown  in  Figure  3- 1(a)  and  the  same  dimensions  as  HPEOS  Experiment  2.  We  then 
simulated  the  response  of  jointed  specimens  having  the  configurations  shown  in  Figures  3- 1(b)  (single 
joint)  and  3- 1(c)  (triple  joint).  The  effect  of  joints  was  established  by  comparing  the  results  of  the 
benchmark  simulation  in  intact  marble  with  the  results  of  simulations  in  jointed  marble. 

Figure  3-13  shows  the  simulated  particle  velocity  histories  at  each  of  the  gage  planes  and  at  the  free 
surface  of  the  intact  marble  experiment  shown  in  Figure  3- 1(a).  The  first  peak  in  the  particle  velocity 
histories  corresponds  to  peak  stress,  and  it  is  approximately  the  same  on  all  the  gage  planes.  This  peak, 
which  also  corresponds  to  the  Chapman-Jouguet  (C-J)  pressure  of  the  explosives,  is  followed  by  a  gentle 
unloading  wave  that  causes  a  gradual  decrease  in  the  particle  velocity  profiles.  This  is  the  Taylor 
unloading  wave,  which  originates  in  the  explosive  gases.  The  second  plateau  in  the  particle  velocity 
records  is  associated  with  the  stress-free  state  in  the  specimen.  It  is  carried  through  the  specimen  by  a 
backward-traveling  release  wave  that  originates  at  the  back  surface  of  the  specimen. 

To  gain  insight  into  the  effect  of  joints  on  wave  propagation,  the  two  experimental  configurations  with 
normally  oriented  artificial  joints,  shown  in  Figure  3- 1(b)  and  3- 1(c),  were  also  simulated.  The  resulting 
particle  velocity  histories  at  the  first  through  the  fourth  gage  planes  are  compared  with  the  results  of  the 
intact  marble  simulation  in  Figures  3-14  through  3-17.  Two  features  associated  with  the  effect  of  joints 
on  the  simulated  particle  velocity  profiles  are  readily  apparent.  The  first  of  these  features  is  depicted  in 
Figure  3-14.  The  particle  velocity  profiles  on  the  first  gage  plane  (ahead  of  the  joint)  are  modified  by 
spikes — one  spike  for  the  single  joint  configuration,  and  three  spikes  for  the  triple  joint  configuration. 
These  spikes,  and  the  corresponding  decrease  in  stress,  result  from  wave  reflections  at  the  rock/joint 
interfaces.  Interestingly,  these  spikes  are  not  quite  as  apparent  on  any  of  the  gage  planes  located  behind 
the  joint  (Figures  3-15  through  3-17),  although  on  the  second  gage  plane  the  velocity  history  associated 
with  the  triple  joint  configuration  shows  characteristics  similar  to  those  observed  on  the  first  gage  plane. 
The  second  feature  associated  with  the  effect  of  joints  on  the  simulated  particle  velocity  profiles  is  the 
delayed  time-of-arrival  (TOA)  of  the  wave  at  the  second,  third,  and  fourth  gage  planes  (all  located 
behind  the  joint).  This  delayed  TOA  is  barely  perceptible  for  the  single  joint  configuration  but  more 
pronounced  for  the  triple  joint  configuration.  The  TOA  delay  can  be  related  to  the  number  of 
reverberations  required  to  bring  the  joint  to  pressure  and  particle  velocity  equilibrium  with  the  test 
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Hgure  3-1 4.  Simulated  particle  velocity  histories  at  the  first 
gage  plane. 


Figure  3-15.  Simulated  particle  velocity  histories  at  the  second 
gage  plane. 
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■igure  3-1 6.  Simulated  particle  velocity  histories  at  the  third 
gage  plane. 


Figure  3-17.  Simulated  particle  velocity  histories  at  the  fourth 
gage  plane. 
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medium.  It  is  not  surprising,  therefore,  that  the  TOA  delay  associated  with  the  triple  joint  is 
approximately  three  times  the  TOA  delay  associated  with  the  single  joint. 

The  same  effects  observed  in  the  PUFF  simulations  can  also  be  observed  in  the  experimental  Hata,  as 
illustrated  in  Figures  3-18  and  3-19  for  the  single  joint  and  triple  joint  conrigurations,  respectively.  In 
each  figure,  the  results  from  an  experiment  on  intact  marble  (HPEOS  Experiments  6  and  7)  arc 
superimposed  onto  the  results  of  an  experiment  on  jointed  marble  (HPEOS  Experiment  10  for  the  single 
joint  and  HPEOS  Experiment  12  for  the  triple  joint).  Part  (a)  in  each  figure  corresponds  to  the  gage 
location  in  front  of  the  jointed  section  of  the  specimen,  whereas  part  (b)  corresponds  to  the  gage  location 
behind  the  jointed  section.  The  narrow,  large  amplitude  spikes  at  t  =  8  |j.s  in  Figure  3- 18(a)  and  at 
t  =  1 1  p,s  in  Figure  3-19(a)  are  most  likely  due  to  an  air  gap  at  the  interface  between  the  joint  and  the 
rock.  This  explanation  is  plausible  because  the  spike  is  not  observed  at  any  of  the  gage  planes  behind 
the  joint  and  also  because  the  spike  reaches  the  gage  before  the  expected  reflection  from  the  joint,  an 
indication  that  the  source  of  the  spike  is  located  at  the  interface  between  the  gage  plane  and  the  joint. 

The  two  features  that  are  associated  with  the  effect  of  joints  on  wave  propagation  in  the  simulated 
particle  velocity  histories  can  also  be  observed  in  the  experimental  records.  On  the  first  gage  plane,  the 
particle  velocity  history  is  modified  by  the  appearance  of  spikes — one  spike  for  the  single  joint 
configuration  and  three  spikes  for  the  triple  joint  configuration.  On  the  second  gage  plane,  the  effect  of 
joints  is  observed  primarily  in  the  form  of  a  delayed  TOA  of  the  wave  at  the  gage  plane.  For  the  triple 
joint  configuration,  a  slight  perturbation  to  the  particle  velocity  history  at  the  second  gage  plane  can  also 
be  observed.  The  larger  than  expected  delay  in  the  TOA  of  the  wave  at  the  second  gage  plane  for  the 
experiment  with  a  single  joint  is  due  to  the  presence  of  an  air  gap  at  the  rock-joint  interface  and,  to  a 
lesser  extent,  to  slight  geometrical  differences  between  the  intact  and  jointed  samples. 

The  results  presented  so  far  pertain  to  the  effect  of  sand-filled  normal  joints  with  rough  surfaces.  Similar 
effects  were  observed  for  sand-filled  normal  joints  with  smooth  surfaces.  A  comparison  of  the  results  of 
HPEOS  Experiment  13  (smooth  triple  joint)  with  the  results  of  HPEOS  Experiment  12  (rough  triple 
joint)  reveals  no  detectable  differences  in  the  response  of  smooth  and  rough  triple  joints.  The  same  can 
be  said  of  the  response  of  smooth  and  rough  single  joints,  except  that  the  particle  velocity  gages  on  the 
first  gage  plane  of  HPEOS  Experiment  1 1  (smooth  single  joint)  failed  before  the  arrival  of  the  reflection 
from  the  joint  at  the  gage  plane. 

Joint  inclination  appears  to  have  had  a  significant  effect  on  the  observed  behavior.  Only  inclined  single 
joints  were  investigated.  The  joint  was  oriented  at  an  angle  of  45°  with  respect  to  the  direction  of  wave 
propagation.  Two  experiments,  one  with  a  rough  joint  (Experiment  9)  and  one  with  a  smooth  joint 
(Experiment  14),  were  conducted,  but  a  failure  in  the  data  acquisition  system  resulted  in  no  data  being 
recorded  in  HPEOS  Experiment  9.  The  data  on  the  first  gage  plane  of  HPEOS  Experiment  14  were 
polluted  by  high  amplitude  spikes  believed  to  be  due  to  an  air  gap  at  the  rock-joint  interface.  However, 
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(a)  Particle  velocity  histories  at  the  first  gage  plane. 


(b)  Particle  velocity  histories  at  the  second  gage  plane. 


Figure  3-1 8.  Effect  of  a  single  sand-filled  joint  on  wave 
propagation. 
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(a)  Particle  velocity  histories  at  the  first  gage  plane. 


(b)  Particle  velocity  histories  at  the  second  gage  plane. 


Figure  3-1 9.  Effect  of  a  triple  sand-filled  joint  on  wave 
propagation. 
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the  particle  velocity  histories  on  this  gage  plane  indicate  that  the  inclined  single  joint  has  the  same 
signature  as  the  normal  joint,  but  with  a  prolonged  duration.  The  signature  of  the  inclined  joint  at  the 
second  gage  plane  resembles  that  of  a  triple  normal  joint  more  than  that  of  a  single  normal  joint.  Both 
these  observations  point  to  the  inclined  joint  as  having  more  effect  on  wave  propagation  than  a  normal 
joint,  all  other  parameters  being  equal.  This  is  not  surprising.  From  Figure  3-1,  it  can  be  seen  that  a 
plane  wave  propagating  through  the  jointed  section  of  a  specimen  with  a  single  normal  joint,  interacts 
with  the  joint  only  as  it  propagates  through  the  joint  itself,  and  thus  the  effective  thickness  of  the  joint  is 
the  same  as  its  actual  thickness.  On  the  other  hand,  a  single  inclined  joint  interacts  with  the  plane  wave 
throughout  the  time  it  takes  for  the  wave  to  propagate  through  the  whole  jointed  layer;  thereby  making 
the  effective  thickness  of  the  joint  equal  to  the  thickness  of  the  layer,  not  the  thickness  of  the  joint.  In 
the  present  study  the  jointed  layer  was  about  13  mm  thick,  compared  with  a  joint  thickness  of  2  mm. 

The  final  joint  parameter  varied  in  this  study  was  joint  thickness.  HPEOS  Experiment  8  was  conducted 
with  a  single  smooth  inclined  joint.  The  joint  was  closed,  i.e.,  the  two  joint  surfaces  were  in  contact. 
The  results  of  this  experiment  do  not  differ  from  the  results  of  the  experiments  on  intact  marble,  an 
indication  that  closed  joints  have  no  effect  on  compressive  wave  propagation. 
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SECTION  4 


SPHERICAL  WAVE  EXPERIMENT 

4.1  MOTIVATIONS  AND  OBJECTIVES. 

Over  the  past  several  years,  DNA  has  sponsored  experimental  programs  that  provided  high  stress, 
dynamic  material  property  data  for  Danby  marble.  These  programs  include  the  large  scale  DISTANT 
MOUNTAIN  test  series  that  provided  data  for  Danby  marble  up  to  stresses  of  about  100  kbar  (Keough 
et  al.,  1993).  The  in-contact  explosive  experiments  described  in  Section  3  of  this  report  also  provide 
data  for  Danby  marble  at  stresses  in  excess  of  100  kbar,  but  at  a  smaller  scale  than  that  of  the  DISTANT 
MOUNTAIN  experiments.  Still  more  experiments  on  Danby  marble  were  performed  by  Furnish  (1994) 
and  Gaffney  and  Smith  (1994),  and  numerous  other  investigations  were  found  in  the  literature.  This  data 
bank  provides  a  wealth  of  information  about  the  behavior  of  marble  under  dynamic  loading,  including 
precursor  and  precursor  decay,  attenuation,  unloading,  and  phase  transitions.  However,  all  these  data  . 
were  obtained  under  one-dimensional  strain  loading  conditions.  The  spherical  wave  experiment 
described  in  this  section  supplements  the  existing  data  by  providing  results  for  marble  in  a  comparable 
stress  range  (~1(X)  kbar)  but  under  spherical  loading  conditions.  Furthermore,  the  experimental  results 
can  be  used  to  validate  and/or  calibrate  constitutive  models  for  use  in  finite  element  and  finite  difference 
hydrocodes. 

4.2  EXPERIMENTAL  TECHNIQUES. 

In  the  spherical  wave  experiment,  a  small  explosive  charge  is  detonated  at  the  center  of  a  27-cm- 
diameter,  27-cm-long  cylindrical  block  of  marble.  The  motion  of  embedded  wire  gages  through  an 
externally  applied  magnetic  field  measures  the  velocity  histories  at  several  radii  from  the  charge 
(Florence  et  al.,  1993;  Gefken  and  Florence,  1993).  The  specimen  was  prepared  from  a  core  obtained 
from  the  same  stock  used  in  the  DISTANT  MOUNTAIN  test  series  as  well  as  for  the  uniaxial  strain 
experiments  of  the  present  study. 

The  overall  configuration  of  the  experiment  is  shown  in  Figure  4-1.  The  specimen  assembly  was  placed 
inside  a  pressure  vessel  and  surrounded  by  a  solenoid  that  applied  a  constant,  vertically  oriented 
magnetic  field  through  the  particle  velocity  gage  planes.  A  6.9-MPa  static  overburden  pressure  was 
applied  to  the  specimen  during  the  experiment.  Beforehand,  the  specimen  was  treated  with  a  water¬ 
proofing  agent  to  ensure  that  it  remains  dry  during  the  experiment. 
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Figure  4-1.  Overall  configuration  of  the  spherical  wave  experiment  in  Danby  marble. 
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A  6.0-g  EL506D  (Detasheet)  charge,  detonated  using  a  mild  detonating  fuse  (MDF),  was  used  to  provide 
the  loading  in  the  spherical  wave  experiment.  Details  of  the  charge  cavity  are  shown  in  Figure  4-2.  As 
shown,  the  EL506D  charge  was  contained  within  a  Lexan  casing,  and  the  MDF  fuze  within  a  stainless 
steel  tube.  The  charge  was  detonated  at  its  center  to  produce  as  symmetric  a  spherical  wave  as  possible 
in  the  rock  sample. 

The  particle  velocity  gages  operated  on  the  same  principle  as  their  counterpart  in  the  uniaxial  strain 
experiments.  The  motion  of  a  gage  element  of  length  i,  moving  through  a  magnetic  field  of  strength  B, 
generates  an  electromagnetic  force  e,which  can  be  related  to  the  velocity  by  Faraday’s  law, 

"'If  ('‘-D 

In  this  experiment,  eight  gages  were  used  to  record  the  radial  particle  velocity  of  the  rock  sample  at 
several  radii  from  the  charge.  The  particle  velocity  gages  consisted  of  strips  of  copper  wire  mounted 
into  precision-machined  grooves  in  the  specimen.  As  shown  in  Figure  4-3,  two  diametrically  opposite 
gages  were  used  at  each  of  eight  radii  starting  at  the  charge  boundary  and  extending  outward  to  a  radius 
of  5.0  cm.  In  most  applications,  the  active  element  of  each  gage  would  cover  the  full  circumference  of 
one  of  the  concentric  circles  shown  in  the  figure.  However,  in  the  present  study  a  different  approach 
was  taken  to  ensure  that  the  direction  along  which  the  particle  velocity  histories  were  measured 
coincided  with  the  direction  of  wave  propagation  in  the  uniaxial  strain  experiments.  This  arrangement 
allows  for  a  direct  comparison  between  the  results  of  the  uniaxial  strain  experiments  and  the  results  of 
the  spherical  wave  experiment  even  if  the  rock  behaves  in  an  anisotropic  fashion. 

4.3  EXPERIMENTAL  RESULTS  AND  DISCUSSION. 

The  particle  velocity  histories  measured  in  the  spherical  wave  experiment  on  Danby  marble  are  shown  in 
Figures  4-4  through  4-11.  Each  of  these  figures  shows  the  two  particle  velocity  histories  measured  at  a 
given  range.  As  was  shown  in  Figure  4-3,  two  colinear  gage  arrays,  one  on  either  side  of  the  high- 
explosive  charge,  were  used  to  measure  the  particle  velocity  histories.  These  gage  arrays  are  designated 
by  the  letters  ‘a’  and  ‘b’  in  Figures  4-4  through  4-11.  The  particle  velocity  histories  at  the  two  ranges 
closest  to  the  charge  show  characteristics  that  are  different  from  the  remainder  of  the  data.  This 
difference  may  be  related  to  premature  gage  failure  due  to  the  high  stresses  generated  by  the  explosive 
charge,  or  it  may  be  related  to  a  transition  from  a  compaction-dominated  behavior  near  the  charge  cavity 
to  a  crack  propagation-dominated  behavior  further  away  from  the  charge.  The  diameter  of  the  charge 
cavity  increased  by  more  than  25%  (from  an  initial  value  of  2.0  cm  before  the  experiment  to  a  final 
value  of  2.55  cm  after  the  experiment),  indicating  significant  compaction  near  the  charge  cavity. 
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Figure  4-2.  Details  of  the  charge  cavity  and  the  detonating 
fuze  used  in  the  spherical  wave  experiment. 
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Particle  velocity  histories  at  the  1 .0-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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Particle  velocity  histories  at  the  1 .5-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 


VELOCITY  (m/s)  VELOCITY  (m/s) 


Figure  4-6.  Particle  velocity  histories  at  the  2.0-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 


Figure  4-7.  Particle  velocity  histories  at  the  2.5-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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Figure  4-8.  Particle  velocity  histories  at  the  3.0-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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Figure  4-9.  Particle  velocity  histories  at  the  3.5-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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Figure  4-10.  Particle  velocity  histories  at  the  4.0-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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Figure  4-11.  Particle  velocity  histories  at  the  5.0-cm  range  in  the 
spherical  wave  experiment  on  Danby  marble. 
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The  particle  velocity  gages  at  the  2.0-cm  range  captured  the  early  time  response  of  the  rock  but  failed 
to  record  the  subsequent  rebound  signal.  The  late  time  failure  of  the  gages  at  this  range  is  a  strong 
indication  that  the  gages  at  the  1.0-cm  and  the  1.5-cm  ranges  are  likely  to  have  failed  during  the 
experiment,  and  that  the  signals  recorded  by  these  gages  are  not  a  reliable  indication  of  the  actual 
response  of  the  rock. 

At  ranges  of  2.5  cm  or  more,  the  particle  velocity  histories  in  the  specimen  were  recorded  for  the 
duration  of  the  experiment.  At  early  times,  these  histories  are  characterized  by  a  sharp  rise  to  peak 
followed  by  an  outward  motion  that  lasted  for  about  10  ps.  This  portion  of  the  velocity  records  is 
associated  with  the  outward-traveling  wave  that  originates  in  the  explosives.  This  early-time  response  of 
the  rock  is  very  reproducible,  as  indicated  by  the  nearly  identical  responses  of  the  gage  arrays  on  either 
side  of  the  explosive  charge  (see  Figures  4-6  through  4-11). 

The  outward  expansion  phase  was  followed  by  a  rebound  phase.  In  previous  experiments  on  nonporous 
rock  (Miller  and  Florence,  1990)  and  on  porous  CaCDs  rock  (Gefken  and  Florence,  1992),  the  rebound 
signal  caused  an  inward  contraction  in  the  motion  of  the  rock  sample  that  was  manifested  by  a  reversal 
of  the  sign  of  the  particle  velocity  histories.  In  the  present  experiment,  the  particle  velocities  recorded  at 
all  ranges  remained  positive  throughout  the  rebound  phase  of  the  signal.  No  specimen  contraction  was 
observed. 

This  important  difference  between  the  rebound  signals  measured  in  the  present  experiment  and  those 
measured  in  other  investigations  can  be  related  to  cracking  of  the  rock  sample.  In  the  Miller  and 
Florence  (1990)  study,  the  explosive  consisted  of  a  3/8-g  PETN  charge,  as  opposed  to  the  more  powerful 
6-g  Detasheet  charge  of  the  present  investigation.  In  the  Gefken  and  Florence  (1992)  study,  the  same 
6-g  Detasheet  charge  was  used,  but  the  limestone  tested  was  approximately  15%  porous.  Extensive 
cracking  was  not  observed  in  either  of  the  two  earlier  studies.  In  contrast,  the  Danby  marble  specimen 
of  the  present  study  was  severely  cracked,  as  shown  in  Figure  4-12.  Two  distinct  networks  of  cracks  can 
be  seen  in  each  half  of  the  specimen.  The  first  network  consists  of  numerous  cracks  em^ating  from  the 
charge  boundary  and  extending  radially  outward  toward  the  free  surface  of  the  specimen.  These  radial 
cracks  extend  throughout  the  gaged  region  of  the  specimen,  with  a  few  cracks  propagating  all  the  way  to 
the  free  smface.  The  second  network  of  cracks  consists  primarily  of  circumferential  cracks  and  does  not 
appear  to  be  symmetric  with  respect  to  the  center  of  the  explosive  charge.  An  additional  crack  that  does 
not  appear  to  be  associated  with  either  of  the  crack  networks  mentioned  so  far  can  be  observed  spanning 
the  whole  specimen  surface.  This  crack  follows  the  path  of  a  preexisting  in  situ  joint  in  the  rock. 

We  believe  that  the  extensive  cracking  in  the  Danby  marble  specimen  caused  degradation  of  the  material 
properties  and  reduced  the  stiffness  of  the  rock.  For  this  reason,  the  specimen  unloaded  along  a 
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Figure  4-1 2.  Experimentally  observed  crack  patterns  in  the  spherical  wave  experiment  on  Danby  marble  (the  markings  on  the  left 
half  of  the  specimen  are  the  grooves  for  the  particle  velocity  gages  and  their  leads). 


“damaged”  modulus  that  was  substantially  lower  than  the  intact  modulus.  Consequently,  much  of  the 
early  time  deformation  was  not  recovered  during  the  contraction,  or  rebound,  phase. 

The  results  of  the  spherical  wave  experiment  also  show  evidence  of  load-induced  anisotropy,  a  response 
mode  usually  associated  with  cracking  in  brittle  materials.  This  anisotropy  is  manifested  as  a  difference 
in  the  recorded  particle  velocity  histories  along  the  two  different  gage  arrays  of  the  spherical  wave 
experiment.  This  occurred  during  the  rebound  phase  of  the  response  and  can  be  seen  in  the  results 
shown  in  Figures  4-8  through  4-11.  A  careful  comparison  of  the  results  presented  in  these  figures  shows 
a  systematic  difference  between  the  late-time  particle  velocity  histories  recorded  by  gages  in  the  ‘a’ 
array  and  those  recorded  by  gages  in  the  ‘b’  array:  a  single  rebound  signal  with  a  duration  of  about 
15-17  p.s  is  observed  on  the  ‘a’  side  of  the  specimen.  In  contrast,  a  two-spike  reboimd  signal  is  observed 
on  the  ‘b’  side  of  the  specimen;  the  first  spike  has  a  larger  amplitude  and  a  longer  duration  than  the 
second  spike,  but  when  the  two  spikes  on  the  ‘b’  side  of  the  specimen  are  combined,  they  have  the  same 
duration  as  the  single  spike  on  the  ‘a’  side.  We  believe  that  these  systematic  differences  are  closely 
related  to  the  effects  of  cracks  on  wave  propagation.  Specifically,  the  lack  of  symmetry  in  the  cracking 
patterns  caused  the  stress  wave  to  become  nonsymmetric,  and  therefore  gages  located  in  different 
regions  of  the  specimen  recorded  different  motions. 

The  time  of  arrival  of  the  initial  disturbance  at  the  successive  particle  velocity  gages,  combined  with 
the  known  location  of  the  gages,  was  used  to  determine  the  initial  wave  propagation  velocity  in  the 
Danby  marble  specimen.  A  least  squares  linear  fit  of  the  results  indicated  that  the  wave  velocity  was 
6.46  mm/ps.  This  fit  is  compared  with  the  experimental  TOA  data  in  Figure  4-13. 

Radial  particle  displacement  histories  were  obtained  by  temporal  integration  of  the  velocity  histories. 
The  results  are  shown  in  Figures  4-14  through  4-21.  Figure  4-14  indicates  that  the  average  maximum 
displacement  at  a  1.0-cm  range  (i.e.,  the  charge  boundary)  was  2.92  mm.  On  the  other  hand,  the  post¬ 
test  measured  displacement  at  the  cavity  wall  was  2.75  mm.  It  may  be  argued  that  the  reasonably  close 
agreement  between  the  two  independent  measurements  lends  credibility  to  the  particle  velocity  history 
measurement  at  the  charge  boundary.  The  good  agreement  does  not,  however,  explain  the  large 
differences  in  the  peak  particle  velocities  recorded  at  gage  locations  la  and  lb. 

The  displacement  histories  can  be  used  to  determine  the  spatial  deformation  gradient  in  the  specimen, 
which,  with  appropriate  assumptions,  can  be  used  to  compute  the  strain  paths  during  the  test.  However, 
because  of  extensive  cracking,  the  relationship  between  the  radial  motion  (measured)  and  the 
circumferential  motion  of  the  specimen  is  not  obvious.  An  alternative  approach  would  be  to  combine 
the  experimental  data  with  a  physically  based  constitutive  model  to  analyze  the  response  of  the  rock  and 
characterize  the  strain  field  in  the  specimen.  Initial  steps  toward  this  type  of  analysis  were  taken  during 
our  investigation.  The  results  of  this  effort  are  described  in  the  next  section. 
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Figure  4-1 3  Distance-time  of  arrival  data  for  determining  the  wave  velocity  in  the 
spherical  wave  experiment. 
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Figure  4-14.  Displacement  histories  at  the  1 .0-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 


Figure  4-15.  Displacement  histories  at  the  1 .5-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 
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Figure  4-1 6.  Displacement  histories  at  the  2.0-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 


Figure  4-17.  Displacement  histories  at  the  2.5-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 
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Figure  4-18.  Displacement  histories  at  the  3.0-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 


Figure  4-19.  Displacement  histories  at  the  3.5-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 
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Figure  4-20.  Displacement  histories  at  the  4.0-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 


Figure  4-21 .  Displacement  histories  at  the  5.0-cm  range  in  the  spherical 
wave  experiment  in  Danby  marble. 
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4.4  CODE  SIMULATIONS. 


Code  simulations  were  performed  to  aid  in  the  interpretation  of  the  results  of  the  spherical  wave 
experiment.  The  simulations  were  performed  in  DYNA2D  (Whirley  and  Engelmann,  1992)  using  a 
newly  developed,  computationally  oriented  constitutive  model  based  on  the  evolution  of  damage  by 
crack  propagation  (Simons  et  al.,  1995).  This  section  briefly  describes  the  main  features  of  the 
constitutive  model  and  presents  the  simulation  results. 

4.4.1  Description  of  the  Model. 

The  model  used  in  this  study  is  based  on  the  activation  and  subsequent  propagation  of  pre-existing 
penny-shaped  flaws.  In  the  model,  we  assume  that  the  heterogeneity  of  geologic  materials  leads  to  local 
stresses  which,  in  general,  are  different  from  the  externally  applied  stresses.  The  local  stresses  are  given 
by  a  statistical  distribution  centered  about  the  externally  applied  stress.  These  local  stresses,  combined 
with  a  fracture-mechanics-based  stress  intensity  factor  formulation,  are  used  to  monitor  the  activation  of 
cracks.  Once  a  crack  is  activated,  its  propagation  is  governed  by  a  stress-  and  crack-size-dependent 
viscous  growth  law,  and  a  limiting  crack  velocity  equal  to  a  fraction  of  the  shear  wave  velocity. 

Crack  propagation  leads  to  damage  and  stiffness  degradation  of  the  rock.  These  effects  are  accounted 
for  in  the  model  by  tracking  the  cracking  strain  tensor  and  a  scalar  damage  parameter.  The  cracking 
strain  is  a  function  of  the  crack  size  and  the  elastic  properties  of  the  rock;  its  separation  from  the  stress- 
producing  elastic  strain  allows  the  model  to  take  into  account  such  important  phenomena  as  stiffness 
degradation,  softening,  and  volume  dilatancy  (i.e.,  bulking).  The  damage  parameter,  which  represents 
the  overall  damage  on  all  the  crack  planes,  is  a  function  of  the  volume  of  cracked  material. 

It  is  used  in  the  model  as  a  relative  measure  of  material  degradation  and  as  a  failure  criterion.  After 
failure,  the  ability  of  the  material  to  support  tensile  loads  in  the  direction  normal  to  the  crack  plane  is 
diminished. 

The  model  is  a  multiplane  model  that  allows  cracks  to  propagate  on  several  planes  with  different  spatial 
orientations.  This  feature  of  the  model  makes  it  possible  to  account  for  crack-induced  anisotropy,  a 
well-documented  aspect  of  brittle  material  behavior.  Even  when  the  initial  crack  distribution  is  random 
and  the  elastic  response  of  the  material  is  isotropic,  cracks  that  begin  to  propagate  do  so  only  on 
preferential  planes  along  which  the  projection  of  the  applied  stress  tensor  satisfies  the  fracture- 
mechanics-based  crack  activation  criterion.  This  inherent  dependence  of  crack  propagation  on  the 
orientation  of  applied  load,  combined  with  the  multiplane  approach,  induces  a  directional  property  into 
the  model  that,  in  general,  renders  an  initially  isotropic  response  anisotropic.  An  advantage  of  the 
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multiplane  approach  is  that  it  affords  a  realistic  representation  of  jointed  rock  media  where  the  initial 
flaws  are  preferentially  oriented  in  one  or  more  directions. 

4.4.2  Simulation  Results. 

After  choosing  reasonable  values  for  the  material  parameters  in  the  model  on  the  basis  of  published  data 
for  marble  and  similar  geologic  materials,  we  attempted  to  simulate  the  spherical  wave  experiment. 
However,  with  the  cracking  model  alone,  we  were  unable  to  obtain  good  qualitative  agreement  between 
the  measured  and  calculated  particle  velocity  histories  near  the  charge  cavity.  The  calculated  particle 
velocity  histories  had  a  higher  peak  and  a  shorter  duration  than  the  measured  histories,  an  indication  that 
the  simulated  response  of  the  specimen  was  more  elastic  than  the  real  response.  To  obtain  good 
agreement  between  the  simulated  and  measured  results,  we  found  it  necessary  to  use  a  porous 
compaction  model  to  represent  the  material  behavior  in  the  vicinity  of  the  charge  cavity  (within  a 
diameter  of  3.5  cm). 

The  simulated  crack  patterns  in  the  spherical  wave  experiment  are  shown  in  Figure  4-22  along  with  the 
experimentally  observed  crack  patterns.  The  model  reproduced  the  radial  and  circumferential  crack 
networks  observed  in  the  experiment.  The  crack  visible  along  the  path  of  the  in  situ  joint  in  the  rock  was 
not  reproduced  in  the  simulation  because  the  joint  itself  was  not  modeled  (i.e.,  initially  isotropic 
behavior  was  assumed  in  the  calculations). 

The  model  was  not  quantitatively  calibrated  to  reproduce  all  the  details  of  the  experimental  data.  The 
objective  of  the  simulations  was  to  try  to  provide  a  better  understanding  of  the  mechanisms  that 
contributed  to  the  observed  response.  To  this  end,  our  results  indicate  that  close  to  the  charge  cavity  the 
behavior  of  the  rock  is  dominated  by  porous  compaction,  whereas  away  from  the  charge  the  behavior  is 
dominated  by  crack  propagation. 
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APPENDIX  A 


RISE  TIME  OF  IN-MATERIAL  GAGES 


Solids  cannot  sustain  discontinuities  in  pressure  or  particle  velocity.  If  a  discontinuity  is  introduced,  the 
material  produces  waves  to  carry  the  discontinuity  away  so  as  to  leave  the  material  in  such  a  state  that 
both  pressure  and  particle  velocity  are  continuous.*  When  a  gage  package  that  has  finite  thickness  and  a 
shock  impedance  different  from  that  of  the  test  medium  is  embedded  between  two  slabs  of  material,  the 
impedance  mismatch  between  the  gage  package  and  the  sample  causes  discontinuity  in  the  flow  field. 
The  perturbed  flow  gives  rise  to  stress  waves  that  originate  at  the  gage  package/specimen  interfaces  to 
carry  away  the  discontinuity.  Several  stress  wave  reverberations  occur  within  the  gage  package  before 
the  gage  comes  to  pressure  and  particle  velocity  equilibrium  with  the  surrounding  medium. 

The  rise  time  of  the  gage  (i.e.,  the  time  required  to  achieve  pressure  and  particle  velocity  equilibrium) 
varies  depending  on  the  gage  plane  thickness  and  the  mismatch  in  material  properties  between  the  gage 
package  and  the  test  medium.  A  good  estimate  of  the  rise  time  of  the  gage  is  important  especially  to 
help  in  the  interpretation  of  experimental  records  when  time-resolved  details  of  the  wave  structure  are  of 
primary  importance.  This  Appendix  presents  a  method  for  calculating  the  rise  time  of  the  gage. 


A.1  ANALYTICAL  SOLUTION. 

Consider  an  in-material,  Lagrangian  gage  embedded  between  two  blocks  of  a  material  sample.  For  the 
purposes  of  the  present  analysis,  the  sample  blocks  can  be  considered  semi-infinite.  Consistent  with  the 
experiments  conducted  during  this  investigation,  we  assume  that  the  average  impedance  of  the  gage 
package,  Z^,  is  smaller  than  that  of  the  sample,  Z^.  We  further  assume  a  linear  acoustic  response. 

If  a  square  pressure  pulse  of  magnitude  Pq  impinges  on  the  gage  plane,  it  will  be  partially  transmitted 
into  the  gage  package  and  panially  reflected  back  into  the  sample  block,  as  shown  in  Figure  A- 1(a). 

The  transmission  and  reflection  coefficients,  T^g  and  may  be  defined  as  follows: 


T  —  ^  ^ 

^^«"Po"Zg-hZg 


(A-1) 


*Haynes,  D.B.,  "Introduction  to  Stress  Wave  Phenomena,"  Sandia  National  Laboratwy  RepcMt  SLA-73-0801,  August  1973. 
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Figure  A-1 .  Wave  interactions  between  the  test  medium  and  the  gage  package  (compressive  waves  are  shown  positive). 


and 


p  _Po*Pl  Zj.Zp 
Po  -z,+z^ 


(A-2) 


where  Pi  is  the  equilibrium  pressure  after  wave  interaction. 

The  next  interaction  of  interest  occurs  as  the  pressure  pulse  of  magnitude  Pi  impinges  on  the  rightmost 
interface  between  the  gage  package  and  the  sample,  as  shown  in  Figure  A- 1(b).  As  before,  the  wave  is 
partially  transmitted  and  partially  reflected.  The  transmission  and  reflection  coefficients,  T^^  and  Kgs, 
are  given  by  the  following  relations. 


and 


T  P2  2  Zy 
^Ss-Pi-Zs  +  Zg 


Kgs- 


P1-P2 

Pi 


Zs  +  Zg 


(A-3) 

(A-4) 


where  P2  is  the  equilibrium  pressure  after  wave  interaction. 

It  is  important  to  note  the  differences  between  Equations  (A-1)  and  (A-2)  and  Equations  (A-3)  and 
(A-4).  These  differences  stem  from  the  fact  that  Equations  (A-1)  and  (A-2)  pertain  to  shock  waves 
traveling  from  a  high  impedance  material  into  a  low  impedance  material,  while  Equations  (A-3)  and 
(A-4)  pertain  to  shock  waves  traveling  from  a  low  impedance  material  into  a  high  impedance  material. 
These  differences  are  further  illustrated  in  Figure  A-1  (i.e.,  the  reflected  pulse  is  a  release  wave  in  one 
case  and  a  compressive  wave  in  the  other  case). 

Figure  A-2  shows  an  x-t  diagram  and  a  P-u  diagram  for  several  wave  reverberations  within  the  gage 
package.  The  pressure  in  the  gage  package  increases  with  each  reverberation.  After  several  wave 
reflections,  the  pressure  in  the  gage  approaches  the  equilibrium  pressure  (Pq).  These  observations  can  be 
formalized  with  the  aid  of  Equations  (A-1)  through  (A-4)  into  a  general  equation  relating  the  relative 
pressure  in  the  gage,  ,  to  the  impedance  of  the  gage  package  and  test  medium,  Zg  and  Zs,  and  to  the 

number  of  reverberations  within  the  gage  package,  n.  This  equation  takes  the  form 


n-1 

1=0 


(A-5) 


The  power  series  in  Equation  (A-5)  is  a  convergent  series  and  possesses  the  following  unique  solution  as 


(A-6) 


PRESSURE  time 


LAGRANGIAN  DISTANCE 

(a)  x-t  diagram  showing  wave  reverberations  in  gage  plane 


PARTICLE  VELOCITY 

(b)  P-u  diagram  showing  the  gage  plane  equilibrium  process 


Figure  A-2.  \Nave  reverberations  due  to  impedance  mismatch  between  the  test  medium 
and  the  gage  plane. 

A- 4 


and  since 


Tsg  +  —  1 


Equation  (A-6)  simply  reduces  to 


(A-7) 


(A-8) 


as  one  would  expect. 

RISE  TIME  ESTIMATE. 

The  rise  time  of  the  gage  is  directly  proportional  to  the  number  of  stress  wave  reflections  within  the  gage 
package  required  to  bring  the  gage  to  pressure  and  particle  velocity  equilibrium  with  the  surrounding 
medium.  If  n  reflections  are  required  to  achieve  equilibrium,  the  rise  time  can  be  obtained  from  the 
relation 


t  = 


n£ 

C 


(A-9) 


where  £  is  the  thickness  and  C  is  the  average  wave  speed  of  the  gage  package.  The  relative  pressure  in 
the  gage  plane  at  time  t  (or  after  n  reverberations)  is  given  by  Equation  (A-5). 

In  the  rock  equation  of  state  experiments  of  the  present  study,  a  typical  gage  package  is  about 
0.0762  mm-thick  and  consists  mainly  of  ejwxy.  In  the  stress  range  of  interest,  the  average  mechanical 
impedance  of  the  gage  package,  Zg,  is  approximately  750  kbar-|xs/cm  and  the  shock  velocity,  C,  is 
approximately  0.6  cm/p.s.*  The  mechanical  impedance  of  the  rock  sample,  Z5,  is  approximately 
1500  kbar-|is/cm.** 

The  equilibrium  pressures  and  corresponding  rise  times  for  the  first  12  wave  reflections  within  a  typical 
gage  plane  are  summarized  in  Table  A-1,  and  the  results  are  plotted  in  Figure  A-3.  As  these  results 
indicate,  three  reverberations  are  required  to  reach  about  95%  of  the  equilibrium  pressure  in  a  time  span 
of  about  38  ns. 


*Keough,  D.,  "Procedure  for  Fabrication  and  Operation  of  Manganin  Shock  Pressure  Gages, "AFWL-TR-68-57,  Air  Force 
Weapons  Laboratory,  Kirtland  AFB,  NM,  19^. 


**Murri,W.  J.,  D.  E.  Grady,  and  K.  D.  Mahrer,  "Equation  of  State  of  Rocks"  Stanford  Research  Institute  Technical  Report 
(SRI  Project  PYU-1883),  July  1975. 
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Table  A-1.  Computed  relative  pressures  and  equilibrium 
times  for  an  epoxy  gage  plane  embedded  in 
limestone. 


n 

Pn 

TV 

t(ns) 

1 

0.6364 

12.70 

2 

0.8678 

25.40 

3 

0.9519 

38.10 

4 

0.9825 

50.80 

5 

0.9936 

63.50 

6 

0.9977 

76.20 

7 

0.9991 

88.88 

8 

0.9997 

101.63 

9 

0.9999 

114.30 

10 

0.9999 

126.98 

11 

0.9999 

139.73 

12 

0.9999 

152.4 
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APPENDIX  B 


SUMMARY  OF  THE  EXPERIMENTAL  RESULTS  OF  THE  UNIAXIAL  STRAIN,  IN¬ 
CONTACT  EXPLOSIVE  EXPERIMENTS 


This  Appendix  summarizes  the  raw  data  and  reduced  experimental  results  of  all  the  uniaxial  strain,  in¬ 
contact  explosive  experiments  conducted  during  the  present  investigation.  Some  of  the  results  presented 
here  have  been  presented  elsewhere  in  this  report.  These  results  are  included  here  for  the  sake  of 
providing  a  complete  set  of  results  that  encompass  all  the  uniaxial  strain  experiments  conducted  during 
the  course  of  the  study. 

Table  B-1  summarizes  the  experimental  conditions  (explosive  system,  joint  type)  of  all  the  uniaxial 
strain  experiments  and  provides  the  thickness  of  each  rock  slab  used  in  constructing  the  specimens  as 
well  as  the  overall  thicknesses  of  the  specimens.  For  each  of  the  experiments  listed  in  Table  B-1,  the 
experimental  results  are  displayed  graphically  in  the  following  order: 

(1)  Specimen  configuration  and  dimensions  (e.g..  Figure  B-4). 

(2)  Summary  of  the  output  of  all  the  stress  or  particle  velocity  gages  in  a  given 
experiment  (e.g..  Figure  B-5). 

(3)  One  figure  for  each  gage  plane  displaying  the  output  of  the  particle  velocity  gages 
(stress  gage  in  the  case  of  HPEOS  Experiment  1)  as  a  function  of  time  (e.g.. 

Figures  B-6  through  B-10). 

(4)  Particle  velocity  histories  (stress  history  in  the  case  of  HPEOS  Experiment  1)  that 
correspond  to  the  gage  output  described  in  (3)  above  (e.g..  Figures  B-1 1  through 


B-1 


Table  B-1 .  Summary  of  the  configurations  and  dimensions  of  the  uniaxial  strain  experiments. 
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(b)  Top  view  of  the  gage  plane 


Figure  B-1.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  1 
(UTTR  limestone). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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Figure  B-2.  Manganin  gage  record  in  HPEOS  Experiment  1. 
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Figure  B-5.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  2. 
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Figure  B-1 1 .  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  2. 


Figure  B-1 2.  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  2. 
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Figure  B-1 3.  Particle  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  2. 


Figure  B-1 4.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  2. 
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Figure  B-1 5.  Particle  velocity  histories  at  the  fifth  gage  plane  (free 
surface)  of  HPEOS  Experiment  2. 
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(b)  Top  view  of  a  typical  gage  plane 


Figure  B-16.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  4 
(marbie). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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Figure  B-17.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  4. 
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Figure  B-1 8.  Output  of  particle  velocity  gages  at  the  first  gage  plane  of 
HPEOS  Experiment  4. 


Figure  B-1 9.  Output  of  particle  velocity  gages  at  the  second  gage  plane  of 
HPEOS  Experiment  4. 
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Figure  B-20.  Output  of  particle  velocity  gages  at  the  third  gage  plane  of 
HPEOS  Experiment  4. 


Figure  B-21 .  Output  of  particle  velocity  gages  at  the  fourth  gage  plane  of 
HPEOS  Experiment  4. 
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(b)  Top  view  of  a  typical  gage  plane 


Figure  B-22.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  5 
(marble). 


Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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Figure  B-23.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  5. 


GAGE  OUTPUT  (V)  GAGE  OUTPUT  (V) 


1.00 


Figure  B-26.  Output  of  particle  velocity  gages  at  the  third  gage  piane  of 
HPEOS  Experiment  5. 


Figure  B-27.  Output  of  particle  velocity  gages  at  the  fourth  gage  plane  of 
HPEOS  Experiment  5. 
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(a)  Overall  front  view  of  specimen 


(b)  Top  view  of  the  first  gage  plane 

Figure  B-28.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  6 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1 . 
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•e  B-29.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  6. 
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Figure  B-30.  Output  of  prefabricated  particle  velocity  gages  at 
the  first  gage  plane  of  HPEOS  Experiment  6. 


Figure  B-31 .  Output  of  SRI-manufactured  particle  velocity  gages  at 
the  first  gage  plane  of  HPEOS  Experiment  6. 
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Figure  B-34.  Output  of  SRI-manufactured  particle  velocity  gages  at 
the  fourth  gage  plane  of  HPEOS  Experiment  6. 


Figure  B-35.  Comparison  of  the  output  of  prefabricated  and  SRI-manufactured  particle  velocity  gages  in 
HPEOS  Experiment  6  (data  tor  the  first  gage  plane). 
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Figure  B-36.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  6  (prefabricated  gages). 


Figure  B-37.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  6  (SRI-manufactured  gages). 
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Figure  B-38.  Particle  velocity  histories  at  the  second  gage  piane  of 
HPEOS  Experiment  6  (SRI-manufactured  gages). 
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Figure  B-39.  Particie  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  6  (SRI-manufactured  gages). 
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Figure  B-40.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  6  (SRI-manufactured  gages). 
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Figure  B-41 .  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  7 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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Figure  B-42.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  7. 
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Figure  B-45.  Output  of  SRI-manufactured  particle  velocity  gages  at  the 
second  gage  plane  of  HPEOS  Experiment  7. 


Figure  B-46.  Output  of  SRI-manufactured  particle  velocity  gages  at  the 
third  gage  plane  of  HPEOS  Experiment  7. 
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Figure  B-47.  Output  of  SRI-manufactured  particle  velocity  gages  at  the 
fourth  gage  plane  of  HPEOS  Experiment  7. 
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Figure  B-48.  Comparison  of  the  output  of  prefabricated  and  SRI-manufactured  particle  velocity  gages  in  HPEOS  Experiment  7 
(data  for  the  first  gage  plane). 
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Figure  B-49.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  7  (prefabricated  gages). 


Figure  B-50.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  7  (SRI-manufactured  gages). 
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Figure  B-51 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  7  (SRI-manufactured  gages). 


Figure  B-52.  Particie  veiocity  histories  at  the  third  gage  piane  of 
HPEOS  Experiment  7  (SRI-manufactured  gages). 
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Figure  B-53.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  7  (SRI-manufactured  gages). 
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Figure  B-54.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  8 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Tabie  B-1 . 
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Figure  B-55.  Summary  of  the  output  of  all  particle  velocity  gages  in 
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Figure  B-56.  Output  of  particle  velocity  gages  at  the  first  gage  plane  of 
HPEOS  Experiment  8. 


Figure  B-57.  Output  of  particle  velocity  gages  at  the  second  gage  plane  of 
HPEOS  Experiment  8. 
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Figure  B-58.  Output  of  particle  velocity  gages  at  the  third  gage  plane  of 
HPEOS  Experiment  8. 


Figure  B-59.  Output  of  particle  velocity  gages  at  the  fourth  gage  plane  of 
HPEOS  Experiment  8. 
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Figure  B-60.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  8. 


Figure  B-61 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  8. 
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Figure  B-62.  Particle  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  8. 


Figure  B-63.  Particie  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  8. 
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Rough,  Sand-Filled  Joint 


(b)  Top  view  of  a  typical  gage  plane 


Figure  B-64.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  10 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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lummary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experimeni 
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Figure  B-68.  Output  of  particle  velocity  gages  at  the  third  gage  plane  of 
HPtoS  Experiment  10. 


Figure  B-69.  Output  of  particle  velocity  gages  at  the  fourth  qaqe  plane  of 
HPEOS  Experiment  10. 
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Figure  B-70.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  10. 


TIME  (ns) 


Figure  B-71 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  10. 
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Figure  B-72.  Particle  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  10. 


TIME  (ps) 

Figure  B-73.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  10. 
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Figure  B-74.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  1 1 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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igure  B-75.  Summary  of  the  output  of  all  particle  velocity  gages  in 
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Figure  B-76.  Output  of  particle  velocity  gages  at  the  first  gage  piane  of 
HPEOS  Experiment  11. 


Figure  B-77.  Output  of  particie  veiocity  gages  at  the  second  gage  plane  of 
HPEOS  Experiment  1 1 . 
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Figure  B-80.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  1 1 . 
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Figure  B-81 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  1 1 . 
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Figure  B-82.  Particle  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  1 1 . 


TIME  (ps) 

Figure  B-83.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  1 1 . 
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(b)  Top  view  of  a  typical  gage  plane 


Figure  B-84.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  12 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1 . 
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ry  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  12. 
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Figure  B-86.  Output  of  particle  velocity  gages  at  the  first  gage  plane  of 
HPEOS  Experiment  12. 


Figure  B-87.  Output  of  particle  velocity  gages  at  the  second  gage  plane  of 
HPEOS  Experiment  12. 
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Figure  B-90.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  12. 


TIME  (ps) 

Figure  B-91 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  12. 
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Figure  B-92.  Particle  velocity  histories  at  the  third  gage  plane  of 
HPEOS  Experiment  12. 


Figure  B-93.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  12. 
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(b)  Top  view  of  a  typical  gage  plane 


Figure  B-94.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  13 
(marble). 

Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1 . 
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Figure  B-98.  Output  of  the  particle  velocity  gages  at  the  third  gage  plane  of 
HPEOS  Experiment  13. 


Figure  B-99.  Output  of  the  particle  velocity  gages  at  the  fourth  gage  plane  of 
HPEOS  Experiment  13. 
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Figure  B-1 00.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  1 3. 


TIME  (ps) 

Figure  B-1 01 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  13. 
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Figure  B-1 02.  Particle  velocity  histories  at  the  third  qaqe  plane  of 
HPEOS  Experiment  13. 
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FigureB-1 03.  Particle  velocity  histories  at  the  fourth  qaqe  plane  of 
HPEOS  Experiment  13. 
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(b)  Top  view  of  a  typical  gage  plane 


Figure  B-104.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  14 
(marble). 


Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1. 
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Figure  B-108.  Output  of  particle  velocity  gages  at  the  third  gage  plane  of 
HPEOS  Experiment  14. 


Figure  B-1 09.  Output  of  particle  velocity  gages  at  the  fourth  gage  plane  of 
HPEOS  Experiment  14. 
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Figure  B-1 1 0.  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  14. 
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Figure  B-111 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  14. 
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Figure  B-112.  Particle  velocity  histories  at  the  third  gage  plane  of 


Figure  B-1 1 3.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  14. 
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Figure  B-1 14.  Specimen  configuration  and  dimensions  for  HPEOS  Experiment  16 
(marble). 


Note:  The  explosive  system  used  in  the  experiment  and  the  thicknesses 
of  the  rock  layers  can  be  obtained  from  Table  B-1 . 


Figure  B-1 1 5.  Summary  of  the  output  of  all  particle  velocity  gages  in  HPEOS  Experiment  16. 
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Figure  B-120,  Particle  velocity  histories  at  the  first  gage  plane  of 
HPEOS  Experiment  16. 


Figure  B-121 .  Particle  velocity  histories  at  the  second  gage  plane  of 
HPEOS  Experiment  16. 
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Figure  B-1 22.  Particle  velocity  histories  at  the  third  gage  plane  of 


Figure  B-1 23.  Particle  velocity  histories  at  the  fourth  gage  plane  of 
HPEOS  Experiment  16. 
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APPENDIX  C 


MODEL  SOLUTIONS  FOR  THE  INTERACTION  OF  STEADY-STATE  COMPRESSION 

WAVES  WITH  A  BOUNDARY 

This  appendix  contains  a  copy  of  the  technical  paper  entitled  “Model  Solutions  for  the  Interaction  of 
Steady-State  Compression  Waves  with  a  Boundary.”  This  paper  was  written  by  M.  Cowperthwaite  and 
will  be  published  in  the  proceedings  of  the  4*  International  Symposium  on  Behavior  of  Dense  Media 
Under  Dynamic  Pressures,  to  be  held  in  Tours,  France,  5-9  June  1995. 
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MODEL  SOLUTIONS  FOR  THE  INTERACTION  OF  STEADY-STATE  COMPRESSION 

WAVES  WITH  A  BOUNDARY 


M.  Cowperthwaite 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


C.l  ABSTRACT. 

Model  solutions  for  the  flows  produced  by  the  reflection  of  steady-state  compression  waves  from 
boundaries  with  different  shock  impedances  are  presented  to  exemplify  flow  features  that  are 
ignored  when  the  incident  wave  is  treated  as  a  shock  discontinuity. 

C.2  INTRODUCTION. 

Advances  in  measuring  particle  velocity  histories  in  shock  wave  experiments  lead  to  an  interest  in 
flow  features  that  are  associated  with  the  structure  of  the  shock.  This  paper  is  concerned  with  such 
flow  features  produced  when  a  steady-state  compression  wave  (SSCW)  with  a  finite  rise  time  is 
reflected  from  a  boundary.  Model  solutions  for  such  flows,  produced  when  SSCWs  with  different 
shapes  are  reflected  from  boundaries  with  different  shock  impedances,  are  constructed  by  the  method 
of  characteristics.  Particular  solutions  are  presented  to  demonstrate  the  explicit  dependence  of  these 
flows  on  the  shape  of  the  SSCW  and  exemplify  flow  features  that  are  ignored  when  the  compression 
wave  is  treated  as  a  shock  discontinuity. 

C.3  MATERIAL  MODEL  AND  CHARACTERISTIC  EQUATIONS. 

For  tractability  in  constructing  our  model  solutions,  we  assume  that  the  material  supporting  an 
incident  SSCW  has  a  linear  constitutive  relationship,  and  a  constant  propagation  velocity  C.  Then 

the  stress  (a),  particle  velocity  (u),  and  specific  volume  (v)  in  the  SSCW  are  related  by  the  equations, 
0  =  PoCu  (C-1)  0  =  PoC2[^l-;J^j  (C-2) 

where  the  density  p  =  v^,  the  subscript  o  denotes  the  material's  initial  state  and  Uq  =  Oq  =  0. 

We  let  t  and  h  denote  time  and  Lagrange  distance,  and  write  the  flow  equations  expressing  the 
balance  of  mass  and  momentum  as 


d\  3u 

at "  ah 


(C-3) 


at 


(C-4) 


For  our  linear  material.  Equations  (C-2),  (C-3),  and  (C-4)  can  be  combined  to  obtain  its  characteristic 
equationsl^l  as. 


ah  _  ^  ^ 
aa  da 
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(C-5) 
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^  _  r-  ^  3u  Vo  ^ 

dp  ap  ■  ip'c^  »'°"s  c- 


(C-6) 


^  *°  denote  the  forward  facing  and  backward  facing 

Interring  Equations  (C-5)  and  (C-6),  and  introducing  the  Riemann  invariants  r  and  st^l  leads  to  the 
CQUEtlOnS, 


:  =  g  +  R(P) 


u  +  — —  =  2r(P)  along  C+ 
PoC 


(C-7) 


t  =  -^+S(a) 


u  -  — —  =  2s(a)  along  C_ 
PqC 


(C-8) 


It  follows  from  Equations  (C-7)  and  (C-8)  that  the  C+  and  C_  characteristics  in  the  (h-t)  plane  are 
straight  lines,  and  that  u  and  o  are  related  to  r  and  s  by  the  equations. 


(C-9) 


a  =  poC  (r  +  s) 


(C-10) 


men  r  and  s  are  known  functions  of  a  and  p,  the  equations  P  =  R-l  (t-h/C)  a  =  8*1  Tt+h/n 

C.4  THE  REFLECTED  FLOW  PROBLEM. 

®  P^^^ced  when  a  SSCW  propagates  to  the  right  and  impacts  another 
different  shock  impedance.  The  SSCW  has  a  finite  rise  time  (T)  and  the  impact 
rSuf  r  ^Srange  position  h  =  L.  For  convenience,  we  also  assume  that  the  bounding  material 
(  M)  has  a  linear  constitutive  relationship,  and  we  denote  its  properties  by  the  subscript  b.  The 
equations  governing  the  flow  in  BM  can  then  be  obtained  by  writing  a  subscript  b  on  the  parameters. 

Equations  (C-7)  and  (C-8).  Figure  C-1  shows  a  typical  (i/T-' 
^^)  diagram  for  such  a  flow  produced  by  the  reflection  of  a  SSCW  from  a  material  BM  with  a 
J* velocity  Cb>C,  when  the  interface  BB'  is  situated  at 

^0"sjdering  Figure  C-1,  we  note  that  OI2C  represents  the  incident 
^P^^sents  the  steady  state  reflected  wave  (SSRW),  n’2  represents  the  region  of 
QQrw produced  by  the  interaction  of  the  reflected  and  incident  waves,  and  n’TW?epresents  the 
SSCW  tr^pitted  into  BM.  The  parallel  straight  lines  emanating  from  OC  are  the  C+  ^ 

incident  SSCW.  The  parallel  straight  lines  emanating  from  IF  with  a  negative 
slope  are  the  C_  charactenstics  m  the  SSRW  and  those  emanating  from  IF  with  a  positive  sIotc  are 
the  C+  charactenstics  in  the  transmitted  SSCW.  pu»iuvc  Mope  are 

Sftof ssrwfrn^  ^  flow  produced  by  the  reflection  of  our 

I?CP  ^  penetration  region  n’2  occurs  because  OI2C  has  a  finite 

nse  time.  The  second,  that  for  linear  matenals,  the  flows  in  the  steady-state  transmitted  and  reflected 
waves  have  the  same  nse  time  as  OI2C  and  are  determined,  respectively,  by  the  flow  conditions 
along  n  and  21 .  Moreover,  m  the  limit  as  the  rise  time  in  the  incident  SSCW  approaches  zero  and 
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Figure  C-1 .  A  nondimensional  Lagrange  distance-time  diagram  for  the 

flow  produced  by  the  reflection  of  a  steady-state  compression 
wave  from  a  boundary  with  a  higher  shock  impedance  when 
the  material  supporting  the  incident  wave  and  the  bounding 
material  (BM)  have  linear  stress  (a)  -  particle  velocity  (u) 
relationships. 
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Cl  Epproaches  01,  it  is  clear  that  the  penetration  region  1T2  will  disappear,  the  incident  wave  will 
become  a  step  discontinuity  represented  by  01,  and  the  reflected  and  transmitted  waves  will  also 
become  step  discontinuities,  represented  respectively,  by  13  and  IW.  For  these  reasons,  we  are 
particularly  i^rested  in  the  penetration  region  and  how  its  flow  is  governed  by  the  shape  of  the 

After  consider  these  qualitative  flow  features,  it  is  convenient  to  define  our  reflected  flow  problem  as 
that  ot  calculating  the  flow  in  the  penetration  region  when  the  flow  in  the  incident  SSCW  is 
prescnbed.  Becau^  we  will  refer  to  Figure  C-1  when  making  such  calculations  by  the  method  of 
characteristics,  we  first  denve  the  equations  for  the  C+  characteristics  in  OITC  and  II'TW  and  the  C 
charactenstics  in  IFRS.  u  ^ 

C.5  THE  EQUATIONS  FOR  THE  CHARACTERISTICS  IN  THE  (t/T  -  h/CT)  PLANE. 

First,  we  dvide  the  first  equation  in  Equation  (C-7)  by  T  and  set  R(B)/T  =  B  to  obtain  the 
equation  for  the  C+  characteristics  in  OI2C  as 


1- JL  a 

T-ct'^P 


0<p<l 


(C-11) 


Next,  we  use  t®  to  denote  the  time  along  the  interface  BB'  and  write  the  equations  for  the  C* 
charactenstics  m  BM  emanating  finom  E'  as 


^  ^  t®  <  tB  <  t® 

T  ChT  ’  ti^t»<^ 


(C-12) 


Combining  Equation  (C-12)  with  the  equation  obtained  by  setting  t  =  t®,  h  =  L,  and  B  =  Bk  in 
Equation  (C-1 1),  then  gives  the  equation  for  the  C+  characteristics  in  E'TW  as 

t  C/ti\  L  f.  C\„ 

T  =  ^(CTj-^CT(l-Qj+Pb  ,  0<Pb^l  (C-13) 

A  procedure,  similar  to  that  used  to  derive  Equation  (C-13),  gives  the  equation  for  the  C 
charactenstics  emanating  from  IF  as 


Jl  ^ 

CT  CT 


0<a<l 


(C-14) 


These  equations  for  Ae  C+  ^  C_  characteristics  will  be  used  to  construct  solutions  to  our  reflected 
inadent  SSCWs  with  different  shapes.  The  flows  produced  by  a  linear  ramp  wave 
will  be  constructed  before  those  produced  by  a  nonlinear  wave.  ^ 

C.6  INTERACTION  OF  A  RAMP  WAVE  WITH  A  BOUNDARY. 

Before  writing  the  ^nations  for  the  ramp  wave,  it  is  convenient  to  derive  the  equations  for  r  and  s 
^  ^  incident  and  transmitted  SSCW.  Combining  Equations  (C-1),  (C-9)  and 

(C-10)  pves  the  equation,  r-s  =  r+s,  which  shows  that,  s  =  0  and  r  =  u,  in  an  incident  SSCW. 
Similarly,  it  can  be  shown  that ,  Sb  =  0  and  rb  =  Ub,  in  the  transmitted  SSCW. 

C.6.1  Flow  in  the  Ramp  Wave. 

We  use  u  and  a  to  denote  the  peak  particle  velocity  and  peak  stress  in  our  incident  ramp  wave  and 
write  Its  equations  for  u  and  a  as  ’ 
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Il'v 

rs'Ti^'-cj 

Equations  (C-1 1)  and  (C-15)  then  allow  us  to  write  the  equations  for  r  and  s  in  the  ramp  wave  as 

r=up  ,  s  =  0  for  0<P<1  (C-16) 

We  can  now  use  the  equations  relating  u  and  a  to  r  and  s,  and  the  boundary  conditions  along  H'  at 
the  interface,  u®  =  Ub  and  a®  =  05,  to  derive  the  equations  governing  the  flows  produced  by  the 
reflection  of  the  ramp  wave  at  the  boundary. 

C.6.2  Conditions  Along  II'. 

Remembering,  that  r  =  uP  in  OI2C  and  that  Sb  =  0  in  ITTW,  Equations  (C-9)  and  (C-10)  allow  us  to 
write  the  boundary  conditions  along  H'  as 

u®=up-s®  =  rb  (C-17) 

o®  =  poC  (up  +  s®)  =  (poC)b  rb  (C-18) 

Combining  Equations  (C-17)  and  (C-18)  gives  the  following  equations  for  s®  and  rb 

where  Z  =  poC/(poC)b. 

We  can  now  calculate  the  flow  in  the  penetration  region  112,  the  SSRW  2rR3  and  the  transmitted 
wave  II'TW. 

C.6.3  The  Flow  in  the  Penetration  Region  II'2. 

To  obtain  the  equations  for  the  flow  in  ir2,  we  first  use  the  condition  that  a  =  P  along  IF  and  set  a  = 
P  in  Equation  (C-19)  to  obtain  the  equation  for  s  in  the  reflected  wave  as 

s  =  u^|^a  ,  0<a<l  (C-21) 

We  then  use  Equations  (C-9),  (C-10),  (C-16),  and  (C-21)  to  obtain  the  equations  for  u  and  o  in  n'2 
as 

and  use  Equation  (C-1 1)  and  (C-1 4)  to  rewrite  these  equations  in  terms  of  h  and  t  as 

luZ  \l  L  UL-h\-| 

“■(l+Z)  [t“ct'^z(ct  JJ 
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(C-25) 


(l+Z)  [CT'^zl^T'CT JJ 

for0^t/T-L/CT<  1  and  0  <  (l^h)/CT  <  1/2. 

C.6.4  The  Flows  in  the  IITW  and  2I'R3. 

The  condition  that  Ub  =  rt  in  a  steady-state  transmitted  wave,  together  with  Equations  (C-20), 
(C-1),  and  (C-13)  give  the  equations  for  ub  and  Ob  in  II’TW  as, 

^  "(l+Z)  LT“Cb|^CT J-CT^^-Cb JJ 


(C-26) 

(C-27) 


The  condition  that  P  =  1  in  a  SSRW,  together  with  Equations  (C-22),  (C-23),  and  (C-14)  give  the 
equations  for  u  and  a  in  2rR3  as 


i  =  [**a^(T*A‘-§)]  (C-29) 

C.6.5  Flows  Produced  by  the  Reflection  of  the  Ramp  Wave  from  a  Boundary. 

We  first  consider  the  steady-state  transmitted  and  reflected  waves  produced  by  the  reflection 
of  our  ramp  wave  from  the  boundary.  Subj^ting  Equations  (C-26)  and  (C-27)  and  Equations 
(C-28)  and  (C-29]^  respectively,  to  tljp  conditions  that  Pb  =  1  ttnd  a  =  1  give  the  equations  for  the 
particle  velocity  (ub)  and  the  stress  (Ob)  along  IT  and  I'R  as 


A  A 

Ub.  .,0b  2Z 

A  -  A  “1+7 

u  a 


(C-30) 


It  follows  from  these  equations,  that  ub  >  u  and  Ob  <  a  when  Z  >  1,  but  that  ub  <  u  and  Ob  >  a  when 
Z  <  1.  Consequently,  the  impact  of  the  ramp  wave  produces  a  reflected  rarefaction  wave  when  poC 
>  (PoC)b,  but  produces  a  reflected  compression  wave  when  poC  <  (PoC)b.  Clearly,  the  equalities  in 
Equation  (C-30)  can  also  be  obtained  from  an  impedance  mismatch  calculation  for  such  linear 
materials  in  the  (u-a)  plane,  when  the  incident,  transmitted  and  reflected  waves  are  treated  as  step 
discontinuities. 


We  now  consider  the  flow  in  the  penetration  region  W2.  For  a  SSCW  with  a  finite  rise  time  T  as 
shown  in  Figure  C-1,  only  the  particles  lying  in  the  range,  I7CT  -  1/2  ^  h/CT  <  IVCT,  enter  the 
penetration  region  along  12,  pass  through  it,  and  then  leave  along  21’.  Across  12  and  2r  the  flow 
derivatives  of  these  particles  have  jump  discontinuities,  but  within  n'2  they  are  constant.  For  this 
case,  the  values  for  the  flow  derivatives  of  u  and  a  in  ir2  follow  froin  Equations  (C-24)  and 
(C-25)  as,  ^(u^)/^(l/T)  =  2Z/(1+Z),  a(u/u)/a(h/CT)  =  -2/(l+Z),  a(a/o)/a(t/T)  =  2/(l+Z),  and 
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9(a/&)/9(h/CT)  =  -2Z/(1+Z).  It  is  clear  from  our  derivation  of  ^nations  (C-24)  and  (C-25)  that 
these  flow  derivatives  are  constant  because  the  incident  SSC!W  is  a  ramp  wave. 

We  next  consider  the  curves  of  constant  u  and  a  in  II'2.  Differentiating  Equations  (C-24)  and  (C-25) 
gives  the  equations  for  the  slopes  of  these  u  and  a  contours,  respectively,  as 


(C-31) 


(C-32) 


It  follows  from  Equations  (C-31)  and  (C-32)  that  the  u  and  a  contours  in  ir2  are  parallel  straight 
lines,  and  that  the  sum  of  the  angles  they  m^e  with  the  (h/CT)  direction  is  7t/4.  When  Z  <  1,  the  u 
contours  in  112  are  steeper  than  21'  but  Ae  a  contours  are  not,  and  consequently  u  decreases  along  21’ 
as  a  increases.  Thus,  while  112  is  covered  by  the  u  contours  that  emanate  from  12,  it  is  covered  by 
the  a  contours  that  emanate  from  12  and  21'.  Alternatively,  when  Z  >  1,  the  a  contours  are  steeper 
than  2r  but  the  u  contours  are  not,  and  a  decreases  along  21  as  u  increases.  112  is  then  covered  by 
the  a  contours  emanating  from  12  and  by  the  u  contours  emanating  from  12  and  21'. 

We  continue  this  treatment  of  our  reflected  flow  problem  by  considering  the  two  limiting  cases  for 
the  boundary,  namely,  a  rigid  wall  and  a  free  surface.  However,  because  of  space  limitations,  we 
will  only  consider  the  reflection  of  a  ramp  wave  from  a  rigid  wall  and  that  of  a  nonlinear  wave  from 
a  free  surface. 

C.7  REFLECTION  OF  A  RAMP  WAVE  FROM  A  RIGID  WALL. 

C.7.1  Conditions  Along  II'. 

Conditions  at  a  rigid  wall  are  obtained  by  letting  Z  — >  0  and  Cb  — >  «>  in  pertinent  equations 
presented  here  earlier.  The  boundary  conditions  for  u  and  a  along  BB'  are  obtajped  as,  u®  =  0  and 
a®  =  2oP,  by  combining  Equations  (C-17)  and  (C-18)  with  the  equations,  s®  =  uP  and  rb  =  0, 
obtained  respectively  by  setting  Z  =  0  in  ^nations  (C-19)  and  (C-20).  The  condition  that  t  — » t®  as 
Cb  — >  oo  from  Equation  (7),  shows  that  the  C+  characteristics  emanating  from  IF  become  parallel  to 
the  h/CT  axis  when  BM  is  treated  as  a  rigid  wall. 

C.7.2  The  Reflected  Flow. 

The  equations  for  u  and  a  in  the  penetration  region  n'2  are  obtained  by  letting  Z  — >  0  in 
Equations  (C-24)  and  (C-25)  as 

|  =  (C-33)  (C-34) 

Consequently,  u  is  constant  along  particle  paths  and  a  is  constant  along  isochrones,  and  as  shown  by 
Equations  (C-31)  and  (C-32),  the  slopes  of  the  u  and  a  contours  in  the  penetration  regions  are, 
respectively,  infinite  and  zero  when  Ae  boundary  is  a  rigid  wall.' 

Setting  P  =  1  and  Z  =  0  in  Equations  (C-22)  and  (C-23)  uives  the  following  equations  for 
2r]^3,  u  =  u  (1-a)  and  a  =  a  (1+a),  which  show  that  u/Tdecreases  from  one  to  zero  while 
a/a  increases  from  zero  to  two  as  a  increases  from  zero  to  one  in  the  SSCW  reflected  from 
a  rigid  wall. 
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C.8  THE  REFLECTION  OF  A  NONLINEAR  WAVE  FROM  A  FREE  SURFACE. 

C.8.1  Flow  in  an  Incident  Nonlinear  Wave. 

Ctescription  Of  an  incident  SSCW.  we  assume  that  it  has  a  point  of 
inflection,  and  for  the  sake  of  tractability  prescribe  its  equations  for  u  and  o  as  ^ 


«  =  JLf  h  f 

^  ^|^t-ctJ-2[t-ctJ 


(C-35) 


^e  temporal  and  spatial  derivatives  of  (u/u)  and  (a/a)  in  this  wave  are  zero  when  B  =  0  and  1  its 
point  of  inflection  is  located  at  P  =  1/2,  and  its  equations  for  r  and  s  are 


r  =  u(3p2-2p3),  s  =  0  for  0^p<l 
C.8.2  Conditions  Along  II’. 


(C-36) 


, .  . — •  — wuiiuiuuiia  lui  u  aiiu  o  along  an  are  obtained  as  o®  =  0  anH  n  -  n  ft  i 
combining  Equations  (C-17)  and  (C-18)  with  the  equations,  sB  =  -S  a?d  A  =  0 '  dSn^' 

charactenstics  emanating  from  IT  collapse  onto  IB'  when  the 

C.8  J  The  Reflected  Flow. 

1!,“  condition  qB  =  0  and  Equation  (C-9)  to  show  that  rB  =  -sB  alona  BB'  and 

then  the  condinon  that  a  =  p  along  IT  to  obtain  the  foUowing  equation  for  s  in  the  leKLw, 

s  =  u  (3  a2  -  2  a3)  (C-37) 

Equations  (C-9),  (C-10),  (C-36),  and  (C-37)  then  give  the  equations  for  u  and  a  in  n'2  as 


?  =  3(p2  +  a2)-2(p3  +  a3) 
?=3(p2-a2)-2(p3-a3) 


(C-38) 


(C-39) 


for  0  <  p  <  1  and  0  <  a  <  1.  Setting  p  -  1  in  these  equations  gives  the  following  equations  for  u  and 
a  in  the  2rR3, 


A  =  1  +a2(3-2a) 
u 


(C-40)  x=l-a2(3-2a) 

a 


(C-41) 


which  show  that  Wu  int^ases  from  one  to  two  while  o/a  decreases  from  one  to  zero  as  a  increases 
from  rero  to  one  m  the  SS  rarefaction  wave  reflected  from  a  free  surface.  CotSiSg  l^uaS^ 

eniiarirt^  f  (C-38)  and  (C-39),  and  with  Equations  (C-40)  and  (C-41)  give  the 

^nations  for  u  and  a  in  ir2  and  2rR3  in  teims  of  h  and  t  but  welill  not  presem  S  equations 
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Because  of  an  interest  in  Lagrange  gage  histories  recorded  in  shock  wave  experiments,  we  will  now 
consider  the  flow  produced  by  the  reflection  of  our  nonlinear  wave  from  a  free  surface  in  more 
detail.  The  equations  obtained  by  differentiating  Equations  (C-35X  (C-38),  and  (C-39)  with  respect 
to  t/T  and  h/CT  show,  in  contrast  to  Ae  flow  produced  by  the  reflection  of  a  ramp  wave,  that  the 
temporal  and  spatial  derivatives  of  um  and  o/ct  are  continuous  across  12  and  21  and  depend  on  both 
t/T  and  h/CT  in  ir2.  Other  significant  features  of  the  flow  can  be  obtained  by  considering  a  set  of 
Lagrange  particle  velocity  histories  and  the  shape  of  the  u  and  a  contours  in  112. 

Such  a  set  of  u/u  versus  t/T  records,  supposedly  obtained  from  four  Lagrange  gages,  Gl,  G2,  G3, 
and  G4,  situated  as  shown  in  Figure  C-3  at  h/CT  =  0.8, 1.0, 1.3,  and  1.5,  is  presented  in  Figure  C-2. 

The  records  from  Gl  and  G2  show  clearly  that  particles  at  Lagrange  positions  h  ^  L/CT  -  1/2  have 
the  same  acceleration  in  both  the  incident  SSCW  and  the  reflected  SS  rarefaction  wave.  The  records 
G2,  G3,  and  G4  show  that  u/u  =  1  when  t/T  =  2  and  thus  indicate  that  this  condition  will  be  satisEed 
by  all  the  particles  entering  n'2.  This  is  verified  in  our  following  treatment  of  the  u  and  a  contours 
in  the  penetration  region. 

The  derivatives,  [9  (t/T)/9  (h/CT)]u  and  [9  (t/T)^  (h/CT)]<y,  are  convenient  for  deriving  properties  of 
the  u  and  a  contours  in  n'2.  Differentiating  Equations  (C-38)  and  (C-39)  gives  the  equations  for 
these  derivatives  in  terms  of  a  and  p  as 


P(l-P)-a(l-«) 

(C-42) 

[a  (h/CT)  J,- 

P  (1-P)  +  a  (1-a) 

(3  («T)  ^  . 

P(l-p)  +  a(l-a) 

(C-43) 

■p  (1-P) -a  (1-a) 

It  follows  from  Equation  (C-42),  that  [9  (l/T)/9  (h/CT)]u  =  1, 0,  -1,  and  0,  respectively,  along  12 
where  a  =  0,  along  IF  where  a  =  P,  along  21'  where  P  =  1,  and  along  22'  where  a  +  P  =  1  and  t/T  =  2. 
Thus  in  ir2, 22'  separates  the  u  contours  connecting  points  on  12  to  points  on  12'  from  those 
coiinecting  points  on  21'  to  points  on  2'!'.  The  u  contours  below  22',  that  emanate  from  points  on  12 
satisfying  the  condition  0  <  P  <  1,  initially  have  the  same  slope  as  these  C+  characteristics  in  OI2C 
and  become  perpendicular  to  the  free  surface  along  12'.  The  u  contours  above  22',  that  emanate  from 
points  on  21'  satisfying  the  condition  0  <  a  <  1,  initially  have  the  same  slope  as  these  C_ 
characteristics  in  2rR3  and  become  perpendicular  to  the  free  surface  along  21'.  Such  a  set  of  u 
contours  is  shown  in  Figure  C-3. 

It  follows  from  Equation  (C-43),  that  [9  (i/T)/9  (h/ClOla  =  1,  «>,  -1,  and  respectively,  along 
21  where  a  =  0,  along  n'  where  oc  =  p,  along  2T  where  P  =  1,  and  along  22'  where  a  -H  p  =  1  and  t/T  = 
2.  Thus  in  IVl  where  11'  is  the  a  =  0  isobar,  the  other  isobars  connect  points  on  12  to  points  on  2T. 
The  a  contours  that  emanate  from  points  on  12  satisfying  the  condition  0  <  P  <  1,  have  the  same 
slope  as  these  C+  characteristics  in  OI2C,  are  parallel  to  H'  along  22',  and  at  points  on  21  satisfying 
the  condition  0  <  a  <  1  have  the  same  slope  as  the  C_  ch^acteristics  in  2rR3.  The  fact  that  the 
isobars  in  n'2  are  perpendicular  to  22',  shows  that  the  a/o  histories  at  the  Lagrange  positions  in  the 
penetration  region  have  a  maximum  at  t/T  =  2  except  at  h  =  L/CT.  Such  a  set  of  a  contours  is  also 
shown  in  Figure  3. 

C.9  CONCLUSIONS. 

Model  solutions  for  the  flows  produced  by  the  reflection  of  SSCWs  from  boundaries  with 
different  shock  impedances  were  constructed  and  used  to  exemplify  flow  features  that  are  ignored 
when  an  incident  SSCW  is  treated  as  a  shock  discontinuity.  In  constructing  these  solutions  by  the 
method  of  characteristics,  both  the  material  supporting  the  SSCW  and  the  b^ounding  material  were 
assumed  to  have  linear  (o-u)  relationships,  and  the  incident  wave  was  treated  both  as  a  linear  ramp 
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Figure  C-2.  A  set  of  nondimensional  Lagrange  particle  velocity  (u/u)  -  time  (t/T) 

profiles  recorded  in  the  flow  produced  when  a  steady-state  compression 
wave  propagating  in  a  material  with  a  linear  stress  (o)  -  particle  velocity  (u) 
relationship  is  reflected  from  a  free  surface. 
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Figure  C-3.  A  nondimensional  Lagrange  distance-time  diagram  showing 
the  stress  (a)  and  particle  velocity  (u)  contours  in  an  incident 
nonlinear  steady-state  compression  wave  and  in  the  flow 
produced  by  its  reflection  from  a  free  surface. 
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wave  and  a  nonlinear  wave  with  a  point  of  inflection.  For  such  solutions,  the  steady  state  flows  in 
the  reflected  wave  and  the  transmitted  compression  wave  have  the  same  rise  time  as  the  incident 
wave  and  are  separated  by  a  penetration  region  where  the  flow  is  unsteady.  As  for  shock 
discontinuities,  the  conation  for  the  steady  reflected  flow  to  be,  either  a  compression  or  rarefaction 
wave,  is  that  the  shock  mpedance  of  the  bounding  material,  be  either  greater  or  less  than,  the  shock 
impedance  of  the  material  supporting  the  incident  SSCW.  Rows  in  the  penetration  region  produced 
by  (1)  the  reflection  of  a  ramp  wave  from  a  surface  with  a  variable  shock  impedance  and  from  a  rigid 
wall,  and  (2)  the  reflection  of  the  nonlinear  wave  from  a  free  surface  were  discussed  to  provide  a 
better  understanding  of  the  dependence  of  the  unsteady  flows  produced  by  such  reflections  on  the 
shape  of  the  incident  SSCW  and  the  boundary  conditions  at  the  interface. 

REFERENCE 

1.  R.  Courant  and  K.  O.  Friedrichs,  "Supersonic  Row  and  Shock  Waves,"  Interscience 
Publishers,  Inc.,  New  York,  1948,  Chapter  HI. 


C-13 


DISTRIBUTION  LIST 

DNA-TR-95-47 


DEPARTMENT  OF  DEFENSE 

DEFENSE  INTELLIGENCE  AGENCY 
ATTN:  DT-1 

DEFENSE  NUCLEAR  AGENCY 
2CYATTN:  ISST 

ATTN:  PM  D  LINGER 
ATTN:  PMAS-TV  F  RENSVOLD 
ATTN:  PMTETREMBA 
ATTN:  WE 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
ATTN:  DTIC/OCP 

FIELD  COMMAND  DEFENSE  NUCLEAR  AGENCY 
ATTN:  FCTN  B  HARRIS-WEST 
ATTN:  NVTV 

FIELD  COMMAND  DEFENSE  NUCLEAR  AGENCY 
ATTN:  FCTT-RQ  B  RISTVET 
ATTN:  FCTT-RQ  E  RINEHART 
ATTN:  FCTTDRBALADI 
ATTN:  FCTTJ  HUGHES 
ATTN:  FCTTS  A  MARTINEZ 
ATTN:  FCTTS  J  LEVERETTE 
ATTN:  FCTTS  LT  COL  AULETTA 
ATTN:  FCTTS  DR  REINKE 
ATTN:  FCTTS  P  THOMPSON 

DEPARTMENT  OF  THE  ARMY 

U  S  ARMY  ENGR  WATERWAYS  EXPER  STATION 
ATTN:  E  JACKSON  CEWES-SD-R 
ATTN:  J  ZELASKO  CEWES-SD-R 

DEPARTMENT  OF  THE  AIR  FORCE 

PHILLIPS  LABORATORY 
ATTN:  PL/SUL 

DEPARTMENT  OF  ENERGY 

EG&G,  INC 

ATTN:  DEILERS 

LAWRENCE  LIVERMORE  NATIONAL  LAB 
ATTN:  FHEUZE 
ATTN:  LEWIS  GLENN 
ATTN:  J  RAMBO 
ATTN:  J  WHITE 
ATTN:  WCMOSS 
ATTN:  RWARD 
ATTN:  TECH  LIBRARY 

LOS  ALAMOS  NATIONAL  LABORATORY 
ATTN:  DAVID  KING 
ATTN:  FREDAPP 
ATTN:  TKUNKLE 
ATTN:  TMCKOWN 
ATTN:  J  FRITZ 
ATTN:  C  MORRIS 

2CYATTN:  REPORT  LIBRARY 
ATTN:  J  N  JOHNSON 


ATTN:  THOMAS  DEY 
ATTN:  TOM  WEAVER 

SANDIA  NATIONAL  LABORATORIES 
ATTN:  DIV9321  W  BOYER 
ATTN:  MIKE  FURNISH 

2CYATTN:  TECH  LIB  3141 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

DEFENSE  GROUP,  INC 

ATTN:  ROBERT  POLL 

ENSCO  INC 

ATTN:  P  FISHER 

JAYCOR 

ATTN:  CYRUS  P  KNOWLES 

KAMAN  SCIENCES  CORPORATION 

2CYATTN:  DASIAC 

KTECH  CORP 

ATTN:  E  SMITH 
ATTN:  FRANK  DAVIES 
ATTN:  LLEE 

LOG  ICON  R  &  D  ASSOCIATES 
ATTN:  J  RENICK 

MAXWELL  LABORATORIES  INC 
ATTN:  DR  E  PETERSON 
ATTN:  J  BAKER 
ATTN:  J  MORRIS 
ATTN:  P  COLEMAN 
ATTN:  S  PEYTON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DAN  PATCH 
ATTN:  JACK  KLUMP 
ATTN:  L  SCOTT 
ATTN:  MARTIN  FOGEL 
ATTN:  MR  M  MCKAY 

SRI  INTERNATIONAL 

2CYATTN:  DR  CURRAN 
ATTN:  DR  JIM  GRAN 
ATTN:  MARKGROETHE 
ATTN:  PDECARLI 

2CYATTN:  THANTOUN 

TECH  REPS,  INC 

ATTN:  FMCMULLAN 
ATTN:  RNAEGELI 

TERRA  TEK,  INC 

ATTN:  W  MARTIN 

TITAN  CORPORATION  (THE) 

ATTN:  A  FREDERICKSON 
ATTN:  S  SCHUSTER 

DIRECTORY  OF  OTHER 

MARYLAND  UNIVERSITY  OF 
ATTN:  RICHARD  DICK 


Dist-1 


